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*■  SUMMARY 

In  response  to  Phase  III  of  Contract  F336L5-70-C-I052,  the 
..  necessary  methodology  has  been  developed  to  provide  an  inde- 

pendent analysis  capability  to  Air  Force  Studies  Analysis  (AFSAl 
for  accurately  simulating  the  flight  paths  of  single  main  rotor 
helicopters.  The  equations  developed  under  Phase  1 of  this  con- 
tract have  been  generalized  to  represent  other  single  main  rotor 
helicopters  that  are  currently  operational  and  those  designed 
for  the  near  future. 

This  report  reviews  the  necessary  equations  for  the  determina- 
tion of  power  required  ii  accelerated  flight  and  a method  for 
determining  the  input  coefficients  based  on  the  helicopter's 
physical  parameters  and  power  required  in  one  g flight.  The 
impact  of  the  type  of  rotor  system  on  the  maximum  thrust  which 
can  be  produced  is  examined  in  light  of  the  available  flight 
test  data.  The  concepts  of  energy  maneuverability  for  a heli- 
copter are  examined  using  the  AH-1G  helicopter  as  an  example. 

^ Several  applications  of  these  equations  are  considered  including 

..  terrain  following,  decelerating  turns,  and  low  speed  maneuvering. 
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INTRODUCTION 


Combat  simulations  of  aircraft  require  knowledge  of  the  flight 
path  trajectories  which  the  aircraft  may  follow.  The  flight 
path  trajectory  is  normally  described  in  terms  of  the  Euler 
angles  (.  , v , f)  which  give  the  orientation  (direction)  of  the 
aircraft  in  space  relative  to  some  inertial  reference  frame  and 
the  velocity  vector  of  the  center  of  gravity.  The  flight  path 
trajectories  which  the  aircraft  may  describe  are  dependent  on 
the  performance  and  maneuver  capabilities  of  the  aircraft  which 
must  be  known. 

There  are  two  basic  types  of  accelerated  flight  conditions  which 
must  be  considered:  constant  energy  maneuvers  and  maneuvers 

involving  a loss  of  airspeed  (kinetic  energy),  altitude  (poten- 
tial energy),  or  both.  Constant  energy  maneuvers  result  in 
flight  at  constant  speed  and  altitude.  This  g-capability  is 
power  limited  and  is  primarily  a function  of  excess  power  avail- 
able and  rotor  disc  loading.  The  maximum  normal  acceleration 
which  can  be  attained  with  constant  energy  flights  is  referred 
to  as  the  sustained  g-capability.  For  definition  of  this  capa- 
bility, power  required  must  be  computed  based  on  speed,  altitude, 
and  type  of  maneuver.  Transient  g-capability  refers  to  the 
normal  acceleration  attainable  by  the  rotor.  However,  to  obtain 
this  normal  acce leration,  a loss  of  airspeed  (kinetic  energy;  or 
altitude  (potential  energy),  or  both  occurs.  For  example,  turns 
may  be  made  at  high  transient  g-levels  and  maintained  as  long  as 
altitude  permits.  In  this  sense,  transient  g's  are  unlike  the 
transient  g's  associated  with  fixed  wing  aircraft  performance 
which  may  only  be  sustained  for  seconds.  The  maximum  transient 
g depends  on  the  rotor  blade  loading  coefficient,  i.e.,  thrust 
divided  by  blade  area  and  dynamic  pressure  at  the  blade  tip.  To 
determine  this  capability,  flight  path  deceleration  (rate  of 
loss  of  energy)  must  be  computed  for  the  given  speed,  altitude, 
g-level,  and  type  of  maneuver.  Once  the  transient  and  sustained 
g-capability  are  known  and  the  orientation  from  the  trimmed 
flight  condition  is  known  for  a specified  g-level,  speed,  alti- 
tude, and  type  of  maneuver,  the  flight  path  trajectory  can  then 
be  defined. 

The  concepts  of  kinetic  and  potential  energy  are  used  to  develop 
the  equations  to  compute  flight  path  trajectories.  The  rotor 
speed  degree  of  freedom  is  not  considered  in  the  model  because 
it  would  unduly  complicate  the  representation  of  the  pilot  and 
would  not  significantly  affect  the  accuracy  of  the  computed 
flight  path  trajectories  in  high  speed  flight.  Detailed  expla- 
nation of  developments  of  these  equations  can  be  found  in  Ref- 
erence 1. 
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ENERGY  MANEUVERABILITY 


Concepts 


The  concepts  of  energy  maneuverability  were  originally  applied 
to  fixed  wing  aircraft  to  identify  the  best  areas  of  operation 
during  combat.  Overlays  of  certain  charts  would  indicate  which 
aircraft  had  an  advantage  from  energy  considerations.  These 
concepts,  when  applicable,  were  applied  to  the  helicopter  to 
determine  the  optimal  areas  of  operation.  These  areas  will  be 
identified  during  the  following  discussion. 

A measure  of  the  energy  state  of  a helicopter  at  any  altitude- 
airspeed-RPM-combination  can  be  expressed  as 

E = GW  h + 1/2  mV2  + 1/2  IS2  (11 


The  last  term  of  the  above  equation  is  the  kinetic  energy  of 
the  rotor.  Since  most  helicopters  normally  operate  at  constant 
RPM,  the  rotor  energy  has  been  assumed  constant  for  this  study. 
A convenient  term  in  units  of  feet  can  be  derived  from  equation  1 
by  dividing  by  the  weight,  GW,  and  neglecting  the  rotor  energy 
term.  This  term,  specific  energy  (Es),  expresses  the  energy 
state  of  a helicopter  as  a function  of  altitude  and  airspeed 
as  follows: 


E 

s 


= h 


(21 


Energy  management  involves  control  of  the  rate  of  transfer 
between  energy  levels.  The  rate  of  transfer  between  energy 
levels  can  be  determined  by  taking  the  partial  derivative  with 
respect  to  time  of  equation  2 which  is  expressed  as  follows: 


ifs  = 3h  + V dV  (31 

at  at  g at 


Remembering  that  work  is  accomplished  in  moving  from  one  energy 
level  to  another  and  that  power  is  the  rate  of  doing  work,  then 
the  left  side  of  equation  3 must  be  equal  to  the  excess  power. 
The  excess  power  can  be  determined  by  taking  the  difference  in 
horsepower  available  and  horsepower  required  at  the  altitude- 
airspeed-g  flight  condition.  Therefore,  energy  rate  can  be  ex- 
pressed as 


P 

s 


3E 


s . 


(HP 


AVAIL 


at 


GW 


HP  ) 550 
req 


3h  + V av 
at  g at 


(41 
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From  equation  4 , a value  of  energy  rate  (Ps)  can  be  computed 
for  any  airspeed-g-altitude  condition  for  a given  amount  of 
power  available. 


The  aircraft's  capability  to  change  energy  states  is  limited  by 
the  power  available.  Excess  power  may  be  used  to  increase  alti- 
tude (potential  energy),  to  accelerate  the  vehicle  (kinetic 
energy),  or  to  increase  the  lift  above  the  one  g value  to  change 
the  aircraft's  direction  of  flight.  To  change  the  direction  of 
flight,  a normal  acceleration  is  required.  The  optimal  areas  of 
operation  can  be  identified  from  the  areas  where  change  of  direc- 
tion (normal  acceleration)  is  accomplished  for  minimum  expenditure 
of  energy.  The  turn  radius  for  a specified  load  factor  is  deter- 
mined from  the  following: 


(5) 


Associated  with  minimizing  turn  radius  would 
rate.  Turn  rate  is  computed  by 


g/n2-l 

V 


be  maximizing  turn 


(6) 


The  derivation  of  equation  5 is  based  on  the  relationship 
that  n = 1/cosi.  A discrepancy  between  bank  angle  (i)  data  and 
theory  (n  = l/cos <f)  is  shown  in  Figure  1 for  the  AH-1G  helicopter 
as  reported  in  Reference  2.  As  a result  of  Figure  1,  the 
following  analysis  was  suggested  (Reference  3 ) to  define  the 
relationship  of  normal  load  factor  and  aircraft  attitudes.  Co- 
ordinated flight  requires  only  that  the  ball  of  the  turn  and 
bank  indicator  be  centered.  This  requires  a zero  net  side  force 
in  the  body  axis  system.  A zero  net  side  force  can  occur  in  the 
presence  of  some  sideslip  as  a result  of  aerodynamic  side-force 
characteristics.  If  pitch  attitude,  angle  of  attack,  and  side- 
slip angle  are  considered,  the  analysis  yields: 

n^cosl*  = cos6  cos^t>  + cose  S^K<)  + (7) 


where 


_ tanfrtan  a 
cosC 


The  complete  derivation  of  equation  7 can  be  found  in  Appendix  A. 
A comparison  of  the  test  data  and  data  from  equation  7 is 
presented  in  Figure  2 as  reported  in  Reference  4.  A large 
scatter  band  exists,  but  the  theory  line  now  passes  through  the 
center  of  the  points. 
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Figure  2.  Comparison  of  Theory  and  Flight  Test  Data  for 
Bank  Angle  Versus  Normal  Load  Factor. 
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Energy  Diagram 

The  energy  diagram  is  constructed  from  energy  rate  at  maximum 
power  (PJ  versus  airspeed  and  altitude  for  a specified  value 
of  normal  acceleration  (n).  Lines  of  constant  specific  energy 
CES)  and  energy  rate  ( Ps ) are  presented  for  the  AH-  1G  helicopter 
in  Figures  3,  4,  and  5.  The  units  of  energy  rate  (Ps>  are  feet 
per  second  and  indicate  that  the  aircraft  has  the  necessary 
energy  to  change  its  energy  state.  For  example,  a P~  value  of 
30  feet  per  second  could  result  in  an  acceleration  of  4.2  knots 
per  second  at  80  knots  or  a vertical  rate  of  climb  of  1800  feet 
per  minute.  Normal  acceleration  could  be  increased  until  the 
power  required  equal"  d the  excess  power  reflected  in  the  Ps 
value  (equation  4'.  Therefore,  the  energy  rate  contours  give  a 
measure  of  the  ability  of  the  aircraft  to  change  its  energy 
state,  by  increasing  altitude  or  airspeed,  or  increasing  load 
factor  or  any  combination  of  the  above.  A minimum  time  to  climb 
trajectory  can  be  determined  by  locating  the  points  of  tangency 
between  energy  (Es)  and  energy  rate  (Ps)  contours  as  detailed  in 
Reference  5.  In  Figures  3 and  4 the  contour  defined  by  Ps=0 
identifies  the  steady  state  operating  boundary  of  the  aircraft. 
The  aircraft  cannot  operate  outside  this  contour  without  losing 
energy,  either  in  the  form  of  altitude,  airspeed,  or  both.  On 
the  left  the  boundary  is  the  hover  ceiling.  On  the  top  the 
boundary  is  the  service  ceiling.  On  the  right,  the  boundary  is 
the  power  limited  airspeed. 

The  energy  diagram  can  be  constructed  for  different  values  of 
load  factor.  An  energy  diagram  for  the  AH-1G  helicopter  for  a 
load  factor  of  1.5  is  presented  in  Figure  5.  The  negative 
values  of  energy  rate  (Ps)  indicate  that  the  power  supplied  by 
the  engine  is  insufficient  and  that  either  altitude  or  airsneed 
or  both  must  be  lost  in  exchange  for  the  total  required  Dower. 

A maneuver  in  the  area  of  negative  energy  rates  (Ps)  uses  the 
rotor's  transient  g caoability  which  may  be  limited  by  blade 
stall,  rotor  instabilities,  pilot  comfort  or  vertigo  at  high 
turn  rates,  or  vibration.  The  power  required  for  transient  g 
maneuvers  is  quite  high  and  energy  is  lost  rapidly.  The  rotor's 
capability  to  absorb  power  in  a cyclic-only  maneuver  diminishes 
because  it  approaches  the  autorotative  flow  state  at  high  angles 
of  attack.  Additional  power  may  always  be  absorbed  by  increasing 
collective  pitch  but  rotor  loads  may  increase. 


Maneuver  Diagram 

The  maneuver  diagram  can  be  used  to  identify  the  region  where 
the  aircraft  can  maximize  change  of  direction  for  the  least 
expenditure  of  energy.  The  ability  to  change  direction  is  ex- 
pressed in  terms  of  maximum  turn  rate  and  minimum  turn  radius 
at  a given  energy  rate.  Turn  rate -velocity  diagrams  showing 
lines  of  constant  energy  rate  (Ps),  constant  turn  radius  (R), 
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Figure  5.  Energy  Diagratr  for  AH-1G  Helicopter  at  9000  Pounds 
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and  constant  normal  acceleration  (n)  for  the  AH-lG  helicopter 
for  gross  weights  of  7500  and  9000  pounds  are  presented  in 
Figures  6 and  7.  For  the  helicopter,  turn  rate  is  a maximum  for 
a given  energy  rate  (Ps)  at  low  airspeeds.  Turn  rate  maximums 
for  a given  energy  rate  are  easily  identified  from  Figures  6 and 
7.  Minimum  turn  radius  occurs  at  the  point  of  tangency  between 
the  constant  energy  rate  contours  and  lines  of  constant  radius. 
Figures  6 and  7 indicate  these  values  occur  a.  low  airspeed. 

For  the  helicopter,  operation  in  these  regions  would  always  re- 
quire altitude  (potential  energy!  or  a positive  value  of  excess 
energy  rate  since  kinetic  energy  is  so  low.  The  variation  of 
power  required  with  airspeed  for  a nelicopter  is  such  that  power 
increases  as  airspeed  decreases  fcr  airspeeds  below  about  50 
knots.  A more  important  parameter  than  turn  radius  to  the  heli- 
copter pilot  is  normal  acceleration.  The  region  where  normal 
acceleration  is  maximum  at  a given  energy  rate  is  identified  on 
Figures  6 and  7.  The  intersection  of  the  zero  energy  rate  line 
with  the  lines  of  constant  normal  acceleration  define  the 
steady-state  g boundary.  Operation  outside  of  this  line  will 
require  a corresponding  loss  of  altitude  or  airspeed  or  both 
to  provide  the  necessary  energy  rate  required.  The  impact 
of  normal  acceleration  and  turn  radius  at  a given  energy  rate 
can  be  determined  from  Figures  8 and  9 for  the  AH- 1G  helicopter. 
These  figures  show  the  airspeed  at  which  the  pilot  can  execute  a 
specified  turn  radius  and  still  retain  an  energy  level. 

These  charts  primarily  reflect  the  sustained  g capability  rather 
than  the  transient  g capability  of  the  helicopter.  Simply  by 
increasing  the  installed  power  the  zero  energy  rate  line  would 
move  up  in  Figures  8 and  9 allowing  a higher  normal  acceleration 
versus  speed.  Since  the  helicopter  can  operate  outside  the  zero 
energy  rate  line  by  trading  energy,  either  potential  or  kinetic 
or  both,  the  maximum  rotor  thrust  which  can  be  produced  is 
important.  The  maximum  rotor  thrust  which  can  be  produced  is  a 
function  of  the  type  of  rotor  system  and  is  discussed  in  a later 
section.  This  maximum  rotor  thrust  boundary  is  equivalent  to  the 
stall  and  structural  boundary  on  a V-n  diagram  for  a fixed  wing 
aircraft.  A copy  of  the  computer  programs  used  to  generate  the 
above  charts  can  be  found  in  Appendix  A. 
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DERIVATION  OF  EQUATIONS 

The  necessary  equations  are  developed  below  for  accurately  simu- 
lating the  flight  trajectories  of  single -main -rot or  helicopters. 
Maneuverability  limitations  associated  with  different  types  of 
rotor  systems  are  considered.  A technique  for  simulating  a 
particular  compound  helicopter  is  presented  as  a guide  and 
the  importance  of  the  assumptions  concerning  the  compound  helicop 
ter  are  emphasized. 

Technical  Approach 

Normal  acceleration  is  induced  in  fixed  wing  aircraft  by  applying 
longitudinal  control  or  flap  deflection  to  change  the  wing  lift 
coefficient.  The  equivalent  procedure  for  helicopters  is  to 
change  main  rotor  blade  lift  coefficient.  This  is  usually 
accomplished  by  longitudinal  cyclic  and  collective  pitch  control. 
Normal  acceleration  may  be  induced  by  the  use  of  either  of  these 
controls,  independently  or  together.  Cyclic  control  is  usually 
used  alone  to  produce  normal  acceleration  at  high  speeds  in 
maneuvers  such  as  pullups,  pushovers,  and  banked  turns.  It 
becomes  difficult  to  induce  g at  low  speeds  by  use  of  cyclic 
pitch  alone,  and  so  both  the  collective  and  cyclic  control  are 
used. 

Two  types  of  maneuvers  were  simulated  as  reported  in  Reference  1 
to  determine  if  any  difference  existed  in  power  required  and 
angle  of  attack  caused  by  control  technique.  The  power  required 
for  a symmetrical  pull-up  and  decelerating  turn  were  compared  for 
the  same  speed  and  g level.  The  resulting  power  differences  were 
small  and  indicated  that  either  control  technique  would  be 
satisfactory  for  defining  the  power  required  in  accelerated 
flight.  Description  of  the  method  of  simulation  of  these 
maneuvers  can  be  found  in  Reference  1. 

Power  Equation  - Pure  Helicopter 

A set  of  closed  form  equations  have  been  determined  for  pre- 
dicting power  required  as  a function  of  the  flight  condition, 
normal  acceleration,  and  certain  physical  parameters  of  the 
helicopter.  The  development  of  these  closed  form  equations  is 
presented  below. 

The  power  required  to  overcome  the  flat- plate  drag  is  expressed 
as : 


D = 0 . 5foV2 


(8) 


T 

i 
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The  f in  equation  8 is  the  equivalent  (0^=1)  flat-plate  drag 
area  which  includes  the  contribution  of  the  fuselage,  wing, 
elevator,  fin,  and  rotor  hub.  D is  the  total  drag  of  the 
helicopter.  The  contribution  of  rotor  induced  power  is  repre- 
sented as  follows: 


HPi 


" T5TT 


(10) 


This  induced  power  is  linearly  dependent  on  the  rotor  induced 
velocity.  Wing  induced  drag  is  not  currently  included  but  may 
be  added  easily  if  a large  wing  is  used  on  the  helicopter.  A 
weighting  factor  (K^)  is  used  on  this  term  to  improve  correla- 
tion with  measured  data  in  the  low  speed  range.  This  factor  is 
expressed  as: 

K.  = 1 - (p  - 0.14)2.14 


1 + -3(l  ~ 71$) 

1.0 


p <0.14 

Vi  _ 4.  • V p > 0 • 14 

The  total  velocity  parallel  to  the  thrust  vector,  V 
expressed  as : f 


(11) 


is 


V = 


DV 


p 2r>A  yV2  + 0.866V  2 

P 


(12) 


The  first  term  on  the  right  hand  side  of  equation  12  is  the 
component  of  the  free  stream  velocity,  V,  parallel  to  the  thrust 
vector.  The  second  term  is  the  induced  velocity  and  is  based  on 
momentum  theory. 


V. 

l 


T 

2c A \/v2  + 0.866V  2 
P 


(13) 


The  induced  velocity  is  computed  by  solving  equation  12  for  Vp. 
Since  V is  usually  much  greater  than  Vp,  a close  approximation  to 
Vf , the  second  term  of  equation  12,  is  given  by  setting  Vp  = 0. 
After  Vp  is  determined , then  is  determined  from  equation  13. 

The  rotor  thrust,  T,  used  in  the  above  equations  is  determined 
from 

T = /d2  + (nGW-Lw)2  (14) 

where  L is  the  wing  lift, 
w ° 
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The  wing  lift  is  based  on  the  following  assumptions: 

1.  Wing  has  fixed  incidence. 

2.  Wing  is  positioned  in  the  downwash  of  the  primary 
lifting  rotor. 

The  wing  lift  equation  is  derived  as  follows: 


L = qA  a a 
w w w w 


•a  = cc  +1  - tan 

w fus  w 


-l/Vi\ 

XT} 


Substituting  equation  16  into  equation  15,  and  assuming  Vj/V  to 
be  a small  angle, 


( vi\ 

L — qA  a + a - ■■■  I 

w 1 w wl  fus  W V / 


The  induced  velocity  of  the  main  rotor  is  approximated  by 

nGW  - L 


Dividing  equation  18  by  V and  substituting  into  equation  17 


- = qA  a [c 

W W W l 


nGW  - L 


cc  +1 
I fus  w 


The  main  rotor  induced  velocity  factor  is  defined  as  that  fraction 
of  the  main  rotor  induced  velocity  which  impinges  on  the  wing  in 
the  vertical  plane.  A main  rotor  induced  velocity  factor  of  0.5 
is  assumed  for  this  wing. 


Now  equation  19  becomes 


Lw  = qAwa 


nGW  - L 


•c  + l 
fus  w 


Solving  equation  20  for  1^, 


qAv/aw  /_  , _.  nGW\ 

|i : Vvj  p 
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Equation  21  is  based  on  linear  theory.  From  a least  squares 
curve  fit  technique,  it  was  found  that  better  correlation  with 
G81  (Reference  6)  computed  values  of  wing  lift  was  obtained  if 
the  wing-rotor  interference  term  was  increased  by  17  percent. 
Therefore,  wing  lift  is  calculated  by, 


L = 
w 


^Vw 


(>-t) 


+ 1 


w 


1 .17nGW\ 

5qX  J 


(22) 


Correlation  between  equation  22  and  C81  values  of  wing  lift 
is  shown  in  Figure  10. 


The  power  required  to  overcome  wing  induced  drag  is: 


L2V 

V 

550eb2q 


(23) 


The  blade  profile  power  is: 


HPP 


(vw*) 


bcR 


(24) 


The  ‘o>  h 1 > and  f2  in  equation  24  represents  the  constants  in 
the  drag  expression  for  the  rotor  blade  profile  at  a Mach  number 
of  0.75fR/Vs.  The  term  is  zero  when  a symmetrical  airfoil 
is  used. 

Drag  increases  rapidly  at  high  Mach  number  which  results  in  an 
increase  in  power  required  by  the  rotor.  This  compressibility 
power  is  a function  of  the  Mach  number  of  the  advancing  blade 
tip,  the  airfoil  section  of  the  tip,  and  the  tnrust  coefficient. 
The  Mach  number  of  the  tip  of  the  advancing  blade  is  given  by 
ScR(l+p)/Vs.  if  this  number  is  below  a critical  value,  the 
compressibility  power  is  zero.  The  critical  value  is  determined 
by  evaluating  MCr  at  an  c of  3.5tc/a  for  the  blade  tip  airfoil 
section.  The  3.5tc/a  corresponds  to  the  average  angle  of  attack 
of  the  rotor  and  tc  = 2C-j'/7 , the  blade  loading  coefficient.  If 
the  velocity  of  the  tip  is  greater  than  the  critical  value,  then 
the  amount  above  Mcr  is  determined  by  the  equation: 

/M  = ■ **'  - - M + 0. 75tc  (25) 

s 
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Figure  10.  Wing  Lift  Versus  Normal  Acceleration 
at  17000  Pounds  and  Sea  Level. 
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The  compressibility  power  is  then  a function  of  the  amount  by 
which  the  tip  Mach  number  exceeds  Mcr  as  calculated  by: 


HPc  = AM3  [0.0033  - AM(  0.022-0. 11AM)] 


(26) 


Equation  26  is  an  empirical  expression  based  on  level  flight 
compressibility  correction  for  an  NACA  0012  airfoil. 

The  thrust  coefficient  also  is  used  to  compute  the  stall  power. 

The  value  of  the  thrust  coefficient  at  which  the  power  increases 
due  to  the  partially  stalled  rotor  is  referred  to  as  the  divergent 
thrust  coefficient  ( (t c J d £v ) * If  the  thrust  coefficient  is  less 
than  [tcldivi  then  the  stall  power  is  a function  of  the  amount 
the  thrust  coefficient  exceeds  the  [t^^iv  as  expressed  by: 

HPSTAU.  * [3‘*L0(tc  - |tc|div) 

The  above  expression  is  empirical  and  was  determined  by  computing 
the  best  fit  to  flight  test  stall  power  data  reported  in  Ref erence  4. 

The  last  contribution  to  the  power  equation  comes  from  considera- 
tion of  vertical  velocity.  This  term  can  be  expressed  as: 


3/2 


(27) 


= VvGW 
3W 


(28) 


If  V„  > 0 then  ” = 0.85  and  if  V *•  0 then  r = 0.80. 
v v 

The  sum  of  the  above  components  compose  the  total  power  required 
to  perform  a specified  maneuver.  This  power  may  be  provided  from 
the  engine,  deceleration,  sink  rate,  or  any  combination  of  these 
sources.  The  first  reported  flight  test  data  in  terms  of  power 
required  to  perform  a maneuver  can  be  found  in  Reference  4 for  the 
AH-1G  helicopter.  Correlation  between  power  data  calculated  using 
the  closed  form  equation  and  the  flight  test  data  is  shown  in 
Figure  11  for  the  AH-1G  helicopter. 

Power  Equations  - Compound  Helicopter 

Simulation  of  flight  paths  for  compound  helicopters  is  complicated 
by  the  addition  of  auxiliary  propulsion.  Many  different  main 
rotor  collective  settings  and  auxiliary  propulsion  combinations 
may  be  used  to  obtain  a trim  flight  condition.  Other  primary 
differences  in  flight  simulation  between  pure  and  compound 
helicopters  are  in  pitch  attitude,  use  of  auxiliary  propulsion 
to  accelerate,  and  the  change  in  normal  acceleration  with  angle  of 
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attack.  The  equations  used  below  for  a compound  helicopter  are 
valid  only  for  certain  assumptions.  However,  this  technique  is 
included  as  a guide  for  simulation  of  different  compound 
helicopters. 

Bell  Helicopter  Company  computer  program  C81  (Reference  6)  was 
used  to  obtain  these  data  for  a compound  helicopter.  The  follow- 
ing assumptions  were  imposed  on  the  simulation  of  the  compound 
helicopter  (Reference  1,  Appendix  A) . For  speeds  less  than  100 
knots,  no  auxiliary  thrust  was  allowed  to  trim  the  math  model  in 
level  flight.  At  100  knots,  the  collective  pitr  i was  set  at  11 
degrees  and  the  math  model  trimmed  through  the  use  of  auxiliary 
thrust  as  its  speed  increased.  The  auxiliary  thrust  schedule 
required  by  the  math  model  to  trim  in  level  flight  at  speeds 
greater  than  100  knots  was  fixed.  Specifically,  variations  in 
auxiliary  thrust  from  that  value  required  to  trim  at  a speed  are 
not  allowed. 


In  the  standard  helicopter,  propulsive  force  is  provided  by 
increasing  collective  and  tilting  the  rotor  forward  to  obtain 
greater  forward  speed.  For  the  compound  helicopter,  propulsive 
force  is  provided  by  auxiliary  thrust  at  speeds  greater  than  100 
knots  since  collective  is  not  allowed  to  increase  with  speed. 

Since  an  auxiliary  thrust  control  is  not  available  in  the  current 
Air  Force  model,  the  propulsive  force  (auxiliary  thrust)  of  the 
compound  helicopter  must  be  specified  independent  of  any  pilot 
control.  The  propulsive  power  term  in  the  power  equation  is  a 
function  of  the  equivalent  flat  plate  drag  area  (f).  By  express- 
ing the  auxiliary  thrust  as  an  equivalent  flat  plate  drag  area 
representing  the  propulsive  force  of  the  main  rotor,  an  equiva- 
lent flat  plate  area  was  determined  for  a specified  speed.  From 
the  above  method,  an  equivalent  flat  plate  drag  variation  with 
speed  was  determined  as  presented  in  Figure  12. 

The  total  power  required  for  the  compound  helicopter  is  then  the 
sum  of  the  individual  power  components  as  defined  in  equations 
9,  10,  23,  24,  26,  27,  and  28.  The  f -variation  to  be  used 
in  equation  9 is  shown  in  Figure  12.  The  limitations  of  this 
technique  should  be  realized  and  care  exercised  when  simulating 
a compound  helicopter  which  varies  auxiliary  propulsion  in  a 
different  manner  than  assumed  above. 

The  efficiency  of  the  propeller  (r ) is  assumed  to  be  constant 

and  equal  to  0.85.  P P 

The  wing  effectiveness  in  normal  acceleration  is  considered  below. 
The  normal  force  is  approximated  by 

nGW  = T + L (29) 


i - 
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Differentiating  equation  15  with  respect  to  fuselage  angle  of  attack, 
and  using  equation  16: 


Q | 

— = qA  a 1 + ■*-—  - tan 

d<J£  M w w d°-£  00^  I 


■lWI 


Assuming  a fixed  incidence  wing  and  assuming  Vj/V  is  a small  angle, 


5Vi  = avi  ai 
c?af  “o FT  5a 


From  equation  18  and  29 


Equation  30  now  becomes 


^i_  1 

~&T  ~ 77SV 


= qA< 


a (l. 

w w \ 


3T/3a, 


Differentiating  equation  29  and  substituting  equation  33 
£ = * qAw*w  (>  - ^)1 


T e first  term  of  equation  34  represents  the  increase  in  thrust, 
awd  the  last  term  represents  the  decrement  of  wing  lift  caused  by 
increased  rotor  induced  velocity. 

From  linear  theory,  the  rotor  thrust  change  with  angle  of  attack 
at  a constant  collective  pitch  is 


If  collective  pitch  is  reduced  in  high  g maneuvers,  a greater 
change  in  angle  of  attack  will  be  required  for  the  same  g load 
than  if  the  collective  pitch  is  fixed.  This  is  one  way  the 
division  of  lift  between  the  wing  and  the  rotor  in  maneuvering 
flight  may  be  controlled. 

Substituting  equation  35  into  33  and  using  twice  the  dVi/dT  deriva- 
tive (which  assumes  fully-developed  rotor  wake  on  the  wing),  the 
wing  contribution  becomes 

dL 

asj  - iVJ1  - jjjr)  (36) 
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The  term  in  parentheses  may  be  thought  of  as  the  wing's  "g- 
ef fectiveness"  which  is  zero  if  p = ^a/4,  and  becomes  increasingly 
positive  at  higher  speeds.  Since  a is  typically  about  2'",  if 
- = 0.065,  the  zero  lift  value  of  p is  0.102,  which  corresponds 
to  a speed  of  about  45  knots.  At  this  low  speed  the  small  angle 
assumption  (equation  17)  is  imprecise  enough  that  an  exact 
solution  of  equation  17  demands  the  use  of  tan-i  (V^/V) . The 
g-ef fectiveness  of  the  wing  would  be  0.75  at  180  knots  and  only 
0.5  at  135  knots.  Thus  it  is  apparent  that  a wing  is  effective 
in  maneuvering  flight  only  at  high  forward  speeds. 


Acceleration  and  Deceleration 

Level  acceleration  and  level  deceleration  arc  integral  components 
of  the  maneuver  capability  of  a helicopter.  The  prediction  of 
deceleration  is  important  in  determining  the  amount  of  power 
which  can  be  gained  to  perform  a transient  maneuver. 

The  time  to  accelerate  from  Vf  to  V2  is  determined  by  the  specific 
excess  power,  Ps , and  the  speed,  V,  as  follows: 


t 


(37) 


i 


where  ^ is  an  efficiency  factor.  This  equation  is  obtained 
from  elemental  power  and  force,  F,  relations 


F 


m 


dV 

'St 


P = F-V  = mV 

Since  power  available  has  an  important  effect  on  acceleration 
time,  care  must  be  taken  in  specifying  the  power  applied  versus 
time  when  correlating  with  flight  test  data.  Level  acceleration 
data  for  the  AH-1G  helicopter  can  be  found  in  Reference  7.  A 
comparison  of  predicted  and  measured  data  for  a level  accelera- 
tion from  40  to  100  knots  is  shown  in  Figure  13.  The  horsepower 
applied  versus  time  used  was  taken  from  a representative  time 
history  in  the  above  report.  An  efficiency  of  1.0  was  used. 


During  deceleration,  it  is  assumed  that  all  the  power  required 
to  sustain  level  flight  is  used  to  generate  a decelerative  force. 
In  other  words  the  power  supplied  by  the  engine  to  the  rotor  is 
zero  and  the  rotor  is  autorotat ing . The  time  to  decelerate  from 
V2  to  V]^  is 


299-099-557 


25 


Figure  13.  Acceleration  Performance  for  AH-1G 
Helicopter 
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where  r is  an  efficiency  factor.  An  efficiency  (r)  of  0.8  is 
used  to  allow  the  pilot  a rotor  overspeed  margin  A comparison 
of  predicted  and  measured  data  for  a level  deceleration  from  113 
knots  to  40  knots  for  the  AH-1G  helicopter  is  shown  in  Figure  14. 
The  measured  data  can  be  found  in  Reference  7.  It  should  be 
noted  that  these  data  ire  for  zero  sideslip  and  higher  decelera- 
tion can  actually  be  obtained  by  sideslipping. 

For  simulating  flight  paths,  the  amount  of  excess  power  available 
determines  the  acceleration  as  follows 


rgP 

” S 

~V~ 


(39) 


If  this  horsepower  required  is  greater  than  the  horsepower 
available  ( is  negative)  then  the  above  equation  will  predict 
the  corresponding  deceleration.  Proper  attention  should  be 
given  to  the  efficiency  factor  (r)  of  1 for  acceleration  and 
0.80  for  deceleration.  The  limitation  of  equation  39  is  that 
as  V approaches  zero  V becomes  infinite.  This  is  discussed 
further  on  page  76. 

Alpha  Equation  - Pure  Helicopter 

The  angle  of  attack  is  needed  for  survivability  and  terrain 
following  studies.  A closed  form  expression  for  the  angle  of 
attack  variation  with  velocity  and  g-level  has  been  determined 
by  empirical  techniques  which  use  the  angles  of  attack  computed 
by  C81  as  target  values.  The  equation  is: 


c 


8. 643  2.478( n-1)  + 


10 . 422( 


n-1)] 


A GW 
(-’V)1,6 


1.5 


(40) 


The  first  term  of  equation  40  represents  the  change  in  a caused 
by  normal  acceleration  with  correction  for  altitude,  velocity, 
gross  weight,  and  drag.  The  second  term  represents  the  variation 
in  angle  of  attack  with  velocity  in  one-g  level  flight.  The 
third  term  accounts  for  the  change  in  a for  climb  (positive)  and 
descent  (negative).  The  last  term  represents  the  change  in  a for 
a pilot-induced  acceleration  or  deceleration.  It  should  be 
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emphasized  that  this  term  is  zero  if  the  acceleration  or  decelera- 
tion is  not  pilot  induced.  A flight  condition  where  this  occurs 
is  when  the  horsepower  required  is  greater  than  the  maximum 
horsepower  available.  This  change  in  a results  from  the  differ- 
ence of  applying  the  force  for  acceleration  or  deceleration  at 
the  center  of  gravity  or  at  the  rotor  hub  (pilot  induced). 


The  constants  in  the  angle  of  attack  expression  were  determined 
to  best  fit  the  data  computed  by  the  C81  computer  program  (Ref- 
erence 6)  for  the  AH-1G  helicopter.  An  example  of  the  correla- 
tion between  angle  of  attack  data  calculated  using  the  equation 
and  the  C81  data  is  shown  in  Figure  15.  It  is  recognized  that 
the  coefficients  of  equation  40  are  for  a specific  helicopter, 
but  may  be  used  for  different  helicopters  since  most  helicopters 
have  similar  attitude  variations  with  the  parameters  in  the 
equation.  If  angle  of  attack  data  are  available  from  flight 
test  for  a specified  helicopter,  then  the  coefficients  of 
equation  40  can  be  modified  to  give  a better  fit  while  retain- 
ing the  basic  form  of  the  alpha  equation. 


Alpha  Equation  - Compound  Helicopter 

A closed  form  expression  for  the  angle  of  attack  variation  with 
velocity  and  g level  has  been  determined  by  empirical  techniques 
which  use  the  angles  of  attack  computed  by  C81  (Reference  6)  as 
target  values.  This  equation  applies  only  to  the  compound 
helicopter  specified  in  Reference  1,  Appendix  A,  but  represents 
the  basic  form  of  the  alpha  equation  for  a pure  helicopter  with 
different  coefficients.  This  equation  is  presented  only  as  a 
guide  for  simulating  a compound  helicopter.  The  equation  is 


a = 11.5671  0.8667(n-l)  + 6.8717(n-l) 


/Tgw- 


(t'V) 


T7E 


- 10.8159V 


+ 2.5778  - tan” ^ 


+ 


Act, 


For  speeds  less  than  100  knots,  / c;  from  pilot-induced 
acceleration  or  deceleration  is  asrfollows: 


For  speeds  greater  than  100  knots,  fia*  from  pilot-induced 
acceleration  or  deceleration  is  as  follows: 


(41) 


(42) 


AT 
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Figure  15.  Fuselage  Angle  of  Attack  Versus  True 
Airspeed  for  Symmetrical  Pullups  and 
Pushovers  at  7500  Pounds  and  Sea  Level. 
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The  effect  of  the  change  of  angle  of  attack  on  auxiliary  thrust 
can  be  expressed  as  follows: 


— r— = 0.0046224(n)L*L22 V2  + 1.52446(n)1,132 V + 337.4(n)0,426  (44) 

The  above  equation  is  based  on  data  computed  by  C81  (Reference  6). 

A transition  between  Aa^  and  Acif ^ (equations  42  and  43)  is  provided 

to  avoid  a discontinuity  between  the  Aaf^  for  pure  helicopter  and 

the  Aaf2  for  compound  helicopter.  The  equation  for  transition  is 

as  follows: 

Ac  = [l  - 0.05  (V-135.04)]  [n  + (0.05)(V-!35.04)Aaf  (45) 

r rl  r2 

The  above  equation  is  valid  between  80  and  100  knots  only.  The 
limitations  of  the  above  equation  should  be  realized  and  care 
exercised  in  using  the  equation  for  simulating  a compound 
helicopter  different  from  the  one  in  Reference  1,  Appendix  A. 

However,  the  above  technique  can  be  made  to  simulate  different 
compound  helicopters  if  a small  amount  of  C81  data  which  repre- 
sent the  specified  compound  helicopter  is  available.  The 
technique  is  still  useful  even  though  it  requires  some  dependence 
on  C81  or  flight  test  data. 

Maneuverability  Limits 

The  g-capability  of  a fixed-wing  aircraft  in  maneuvering  flight 
is  usually  determined  by  consideration  of  structural,  aerodynamic, 
and  power  limitations.  The  same  limitations  apply  to  a heli- 
copter, in  a slightly  different  manner.  The  structural  limita- 
tion for  the  helicopter  is  primarily  determined  by  vibration  and 
fatigue  life  of  rotor  system  components.  The  resulting  limit  is 
referred  to  as  the  Vne  speed  (never-exceed  velocity). 

The  aerodynamic  limits  for  the  helicopter  are  not  as  straight- 
forward as  those  for  a fixed-wing  aircraft.  Once  stall  is 
encountered  in  a fixed-wing  aircraft,  lift  no  longer  increases 
with  increasing  angle  of  attack.  Stall  is  not  as  well  defined 
in  a helicopter  as  in  a fixed-wing  aircraft.  Blade  stall  is 
encountered  on  some  parts  of  the  rotor  while  other  parts  will  not 
be  stalled.  Because  stall  has  different  effects  on  different 
rotor  systems,  the  "stall  limit"  can  vary  widely  among  rotor- 
craft.  An  attempt  was  made  to  define  the  "stall  limit"  for  the 
different  rotor  systems  based  on  data  available  from  published 
literature . 


The  data  for  each  rotor  system  are  reduced  to  a rotor  blade 
loading  coefficient  (tc)  versus  advance  ratio  (p)  curve.  This 
coefficient  is  calculated  by  dividing  rotor  thrust  (which  is 
approximately  equal  to  the  load  factor  times  the  weight)  by  the 
blade  area  (bcR)  and  the  average  dynamic  pressure  at  the  blade 
tip  0. 5P  ( AR)2  : 
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c 0.5KrR)2  ( bcR ) 

A brief  discussion  of  the  effect  of  the  type  of  rotor  system  on 
the  thrust  load  limit  follows. 

Systems  that  use  moveable  control  surfaces  on  the  main  rotor 
blade  are  susceptible  to  partial  loss  of  control  and  blade 
torsional  instabilities  when  producing  high  thrust.  A single 
rotor  helicopter  which  uses  this  type  of  system  is  the  HH-2C. 

Based  on  data  from  Reference  8,  blade  stall  may  be  encountered 
if  the  tc-p  curve  indicated  in  Figure  16  is  exceeded.  Intentional 
approaches  to  or  inducement  of  retreating  blade  stall  is  a 
prohibited  maneuver  as  stated  in  the  flight  manual. 

Fully  articulated  rotors  exhibit  large  flapping  and  lead-lag 
blade  motion  and  are  susceptible  to  blade  torsional  instabilities 
(Reference  9)  when  subjected  to  high  thrust.  Several  helicop- 
ters use  this  type  of  rotor  system.  For  the  GH-3C  helicopter, 
blade  stall  may  be  encountered  if  the  tc-t.  curve  shown  in  Figure 
17  is  exceeded  based  on  data  taken  from  Reference  10.  The  flight 
manual  states  that  control  difficulties  will  arise  if  the  stall 
is  allowed  to  fully  develop.  For  the  CH-53  A/D  helicopter,  the 
tc-^  curve  shown  in  Figure  18  represents  the  flight  condition 
where  stall  may  be  encountered  based  on  data  of  Reference  11. 
However,  the  primary  indication  of  blade  stall  is  shown  on  a 
cruise  guide  indicator.  These  values  indicate  the  degree  of 
acceptable  vibration  that  will  sustain  component  service  life. 

The  0H-6A  helicopter  uses  this  type  of  system  also.  The  tc-p 
curve  is  shown  in  Figure  19  for  blade  stall  on  the  0H-6A  helicop- 
ter based  on  data  from  Reference  12.  The  above  data  which  were 
taken  from  respective  flight  manuals  are  assumed  to  be 
representative  data  of  the  particular  type  of  rotor  system. 

However,  it  is  recognized  that  the  data  may  be  conservative. 

The  thrust  load  limit  of  hingeless  rotor  systems  has  not  been 
very  well  defined  because  of  limited  flight  test  experience.  To 
date,  the  limit  that  is  most  often  encountered  is  a structural 
load  limit  in  some  component  of  the  rotor  system. 

Teetering  rotor  systems  may  induce  a vibratory  response  of  the 
fuselage  at  high  thrust.  These  systems  rely  on  a rotor-pylon 
isolation  system  to  minimize  this  vibration  and  show  a high 
tolerance  for  blade  stall.  The  maximum  thrust  achieved  for  the 
AH-1J  helicopter  expressed  in  terms  of  tc  versus  g is  presented 
in  Figure  20  (References  13  and  14).  The  stall  alleviating 
effects  of  pitch  rate  (Reference  15)  minimize  the  oscillatory 
thrust  excitation  especially  for  rotors  with  high  flapping 
inertia.  Thus,  vibration  has  been  shown  by  flight  test  to 
decrease  with  increased  normal  acceleration  at  advance  ratios 
above  0.3  (Reference  16). 
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NOTE:  OPEN  SYMBOLS  - SLS 

FLAGGED  SYMBOLS  - 4000  FT 


BLADE  STALL  LIMI 
0 GW  = 10000 
A GW  = 12000 


AIRSPEED  LIMITS' 
• GW  = 10000  LB 
AGW  = 12000  LB 


DATA  BASIS:  REFERENCE  8,  PAGE  11-13 
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Figure  16.  Rotor  Thrust  Limits  for 
HH-2C  Helicopter 
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Figure  17.  Rotor  Thrust  Limits  for 
CH-3C  Helicopter 
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Figure  18.  Rotor  Thrust  Limits  for 
CH-53A  Helicopter 
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NOTE: 


FLAGGED  SYMBOLS  - 6000  FT 
OPEN  SYMBOLS  - SLS 


BLADE  STALL  LIMITS' 
O GW  = 2600  LB 
O GW  = 2200  LB 


AIRSPEED  LIMITS 

;GW  = 2400  LB  - 
GW  = 2000  IB 
• GW  = 1600  LB 


DATA  BASIS:  REFERENCE  12,  PAGE  14-44 


ADVANCE  RATIO,  \i 


Figure  19.  Rotor  Thrust  Limits  for 
0H-6A  Helicopter 
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The  above  data  on  the  different  types  of  rotor  systems  indicate 
the  average  blade  loading  which  should  not  be  exceeded  for  the 
various  reasons  stated  above.  A tc-p  curve  can  be  used  as  a 
basis  for  comparison  of  the  transient  g-capability  of  different 
rotors  and  helicopters  by  taking  the  ratio  of  maximum  tc,  for 
the  speed  and  type  of  rotor,  to  the  tc  in  one-g  flight  at  the 
altitude  of  interest. 

Simplified  Maximum  Thrust  Model 

The  above  data  indicate  that  the  thrust  limitations  of  different 
rotor  systems  are  based  on  considerations  of  vibration,  control- 
lability, and  fatigue  life.  In  the  event  that  flight  test  data 
are  not  available  to  identify  this  limit,  a method  is  needed  to 
give  a reasonable  approximation.  An  aerodynamically  limited 
value  of  thrust  coefficient  (tc)  at  different  advance  ratios  (p) 
can  be  computed  if  the  variation  of  maximum  lift  coefficient  with 
Mach  number  for  the  blade  section  is  known.  The  mathematical 
model  for  this  computation  is  illustrated  in  Figure  21.  In  the 
reverse  flow  region,  zero  rotor  lift  is  assumed.  In  the  retreat- 
ing blade  region,  the  blade  sections  are  assumed  to  be  at  the 
maximum  lift  coefficiet  t for  the  local  Mach  number.  The  lift  in 
the  advancing  blade  region  is  computed  using  the  maximum  lift 
coefficient  out  to  the  blade  station  at  which  the  resulting  lift 
moments  of  the  advancing  and  retreating  blades  are  equal  and 
opposite.  The  lift  is  assumed  to  be  zero  outboard  of  that 
station.  The  effect  of  pitch  rate  on  the  rotor  is  accounted  for 
by  increasing  the  pitching  moment  which  the  rotor  must  balance. 
This  causes  a reduction  of  stall  on  the  retreating  blade  and  an 
increase  of  thrust  on  the  advancing  blade. 

The  results  of  this  simplified  mathematical  representation  are 
shown  in  Figures  22  and  23  for  the  AH-1J  helicopter  and  the  CH-3C 
helicopter.  The  results  of  this  model  are  not  satisfactory  for 
predicting  the  maximum  thrust  for  a teetering  rotor  system  such 
as  used  by  the  AH-1J.  For  the  articulated  rotor,  a higher  thrust 
was  predicted  than  was  obtained.  Several  explanations  can  be 
offered  for  the  lack  of  correlation.  One  such  explanation  is 
that  two-dimensional  CL-Mach  numbers  were  used.  Another  reason 
is  that  the  articulated  rotor  is  limited  by  other  problems 
before  reaching  its  maximum  thrust  capacity.  Unsteady  aerodynamic 
effects  are  known  to  increase  the  attainable  maximum  lift 
coefficient  of  an  airfoil.  The  conditions  under  which  this 
occurs,  i.e.  , high  rate  of  change  of  angle  of  attack,  do  occur 
on  the  rotor  but  are  not  considered  in  this  simplified  model. 

A derivation  of  the  equations  used  in  the  simplified  math  model 
can  be  found  in  Appendix  B along  with  a copy  of  the  computer 
program. 
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Figure  21.  Analytical  Model  for  Maximum  Rotor  Thrust 
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NOTE:  BASED  ON  VERSUS  MACH  NUMBER 

FOR  SYMMETRICAL  009.3  AIRFOIL 
PITCH  RATE  - 20  DEG/SEC 


NOTE:  DATA  BASIS  FOR  FLAGGED 

SYMBOLS  - REFERENCE  14 


O SYMMETRICAL  PULL-OUT 
□ ROLLING  PULL-OUT 
A CLIMBING  TURN 


ANALYTICAL  MODEL 


DATA  BAS  IS  : REFERENCE  13 
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Figure  22.  Maximum  Thrust  Predicted  by  Simplified 
Analytical  Model  for  AH-1J  Helicopter 
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DETERMINATION  OF  INPUT  COEFFICIENTS 


The  power  equations  were  developed  specifically  for  predicting  the 
total  power  required  for  the  AH-1G  helicopter.  However,  these 
equations  are  capable  of  predicting  the  total  power  required  for 
different  single-main-rotor  helicopters  within  the  accuracy 
required  for  flight  path  simulation.  A set  of  input  coefficients 
can  be  determined  from  the  helicopter's  speed-power  data  at 
different  gross  weights  and  from  certain  physical  parameters 
of  the  helicopter.  Below  is  a discussion  and  example  of  a 
method  to  determine  the  input  coefficients  for  the  power  equations. 

Sensitivity  of  Parameters 

The  limitations  of  the  power  equations  in  representing  compli- 
cated rotor  aerodynamics  for  different  rotor  systems  should  be 
appreciated.  In  seme  cases,  the  equations  may  not  be  capable  of 
representing  the  flight  test  data.  In  general,  the  correlation 
between  flight  test  data  and  computed  data  is  good.  The  follow- 
ing method  is  only  a tool  to  arrive  at  an  initial  starting 
point.  Care  must  be  exercised  in  choice  of  the  input  coefficient 
if  the  values  from  the  INPUT  program  are  not  reasonable.  By 
making  small  changes  in  the  input  data,  correlation  with  flight 
test  data  can  be  improved. 

In  order  to  make  the  necessary  changes  in  the  input  data,  the 
influence  of  each  parameter  on  the  power  curve  must  be  understood. 
A discussion  of  each  of  the  input  parameters  is  given  below.  The 
AH-1G  helicopter  parameters  are  used  in  the  graphic  representa- 
tions . 

1.  Equivalent  Flat-Plate  Drag  (f):  This  parameter  represents 

the  contribution  of  the  fuselage,  wing,  elevator,  fin,  and  rotor 
hub  to  the  parasite  drag  term.  The  change  in  horsepower  with 
speed  and  f is  presented  in  Figure  24.  For  example,  by  increas- 
ing f by  5 square  feet,  the  horsepower  curve  can  be  increased  by 
142  at  140  knots  to  match  flight  test  data  while  increasing  the 
horsepower  at  60  knots  by  only  10. 

2.  Constant  Term  in  Rotor  Blade  Drag  Equation:  C0  is  the  con- 

stant in  the  drag  equation  Cp  = 50  + &2a  • This  term  accounts  for 
the  profile  drag  of  the  rotor  blade  at  zero  angle  of  attack. 
Increases  in  lQ  primarily  result  in  a constant  shift  of  the 

power  curve.  However,  the  effect  of  speed  on  60  is  also 
important  and  is  presented  in  Figure  25.  If  the  calculated 
power  differs  from  the  flight  test  data  by  an  approximately 
constant  value,  then  a change  in  &0  would  improve  correlation. 


3.  Angle  of  Attack  Term  in  Rotor  Blade  Drag  Equation:  62  is 

the  angle  of  attack  coefficient  in  the  drag  equation  C^=tQ+^2a^‘ 
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TRUE  AIRSPEED  - KNOTS 


Figure  25.  Variation  of  Profile  Power  with  True 
Airspeed  and  Blade  Profile  Drag,  60 
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This  term  accounts  for  the  profile  drag  of  the  rotor  blade  with 
angle  of  attack.  The  speed  effect  on  o?  is  presented  in  Figure 
26.  The  difference  in  power  between  different  gross  weights  in 
hover  comes  from  the  contribution  of  62  and  changes  in  induced 
power.  If  the  calculated  power  between  two  different  gross 
weights  differs  from  flight  test  data,  then  a change  in  £>2  w°uld 
improve  correlation. 

4.  Rotor  Induced  Velocity:  The  variation  of  induced  velocity 

with  airspeed  is  presented  in  Figure  27.  A weighting  factor  is 
used  to  increase  the  rotor  induced  velocity  at  low  airspeeds 

to  improve  correlation  with  measured  data.  Usually  this  factor 
is  equal  to  one  above  a speed  of  60  knots.  This  factor  varies 
with  airspeed  and  is  valid  for  several  single  rotor  helicopters. 

5.  Compressibility  Critical  Mach  Number:  Decreases  in  critical 

Mach  number  result  in  increased  power  required  at  high  speeds 

as  shown  in  Figure  28.  Choice  of  the  two-dimensional  critical 
Mach  number  for  a given  blade  section  may  be  a good  initial 
value,  but  the  final  value  may  differ  as  a result  of  not  knowing 
the  lift  coefficient  at  which  the  blade  is  operating  and  not 
knowing  the  rotor's  aerodynamic  environment.  If  the  calculated 
power  does  not  increase  as  rapidly  as  the  flight  test  data  and 
the  tip  Mach  number  is  high,  then  decreasing  the  critical  Mach 
number  will  result  in  improved  correlation. 

6.  Stall:  The  prediction  of  stall  actually  requires  a very 

sophisticated  analysis  which  calculates  the  angle  of  attack  at 
each  radial  and  azimuth  station.  The  method  for  predicting 
stall  in  these  power  equations  is  empirical  and  is  a function  of 
a divergent  value  of  tc  versus  p.  The  amount  of  stall  for  a 

giv  en  flight  condition  depends  on  the  amount  by  which  the  criti- 
cal value  of  tc  is  exceeded.  Therefore,  the  critical  question 
becomes  the  determination  of  this  divergent  tc  versus  p curve 
for  each  helicopter.  From  the  flight  test  data,  stall  is  usually 
seen  as  a rapid  rise  in  power  required  at  high  speed.  If  the 
tip  Mach  number  is  not  high  and  the  flight  test  data  increases 
suddenly  at  high  speeds,  then  lowering  the  divergent  tc  versus 
p curve  would  improve  correlation.  The  increase  of  stall  horse- 
power is  a function  of  the  amount  t„  exceeds  the  divergent  value 
of  tc  as  presented  in  Figure  29.  The  divergent  tc  versus  p 
curve  for  the  AH-1G  helicopter  is 


°div  /l  + 50p2  (47) 

Very  recent  data  have  indicated  that  blade  flapping  inertia 
significantly  affects  the  divergent  tc  versus  p curve.  This 
effect  should  be  investigated  further. 
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Figure  26.  Variation  of  Profile  Power  with  True 
Airspeed  and  Blade  Profile  Drag,  &2 


299-099-557 


46 


BELL-  HELICOPTER  company 


100  120  140 


TRUE  AIRSPEED  - KNOTS 

Figure  28.  Variation  of  Compressibility  Power  with 
True  Airspeed  and  Critical  Mach  Number 
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Methodology 

The  INPUT  program  approximates  the  coefficients  to  the  power 
required  equations  based  on  the  given  power  required  versus  air- 
speed data.  The  correlation  between  predicted  and  actual  values 
of  power  is  given  in  terms  of  percent  error.  The  calculated 
coefficient  may  be  refined  further  by  the  user  to  improve  corre- 
lation. It  is  the  intent  of  this  program  to  only  approximate  these 
values  while  the  final  coefficients  are  determined  by  the  user. 

A flow  diagram  of  the  logic  used  for  this  program  is  presented 
in  Figure  30.  A copy  of  this  program  is  in  Appendix  C.  The 
required  data  for  the  program  are:  (1)  Power  required  versus 
airspeed  for  two  different  gross  weights  at  the  same  altitude, 
density,  and  rotor  rpm;  (2)  hover  power  required  out  of  ground 
effect  for  the  same  conditions  in  (1);  and  (3)  physical  param- 
eters of  the  helicopter:  rotor  radius,  rotor  chord,  number  of 

blades,  rotor  tip  speed,  and  estimated  value  of  equivalent  flat 
plate  drag  area. 

Power  data  at  three  speeds  and  two  weights  are  required  for  the 
actual  computation  of  the  input  coefficients.  These  are  hover, 
minimum  power,  and  maximum  speed  desired.  ^t  each  of  these 
speeds  a residual  power  is  computed.  This  power  results  from 
the  subtraction  of  each  component  of  the  power  equation  as  they 
are  determined  from  the  total  power.  The  remaining  power  repre- 
sents the  contributions  of  the  components  not  yet  determined. 

First,  the  induced  power  is  calculated  and  subtracted  from  the 
total  power  for  both  gross  weights.  Now  the  difference  in  power 
at  hover  between  the  two  gross  weights  is  due  to  the  62  power 
term.  Knowing  this  difference,  t2  is  computed.  The  contribution 
of  £2  t0  power  is  now  subtracted  from  the  residual  power  term. 

Next  the  values  of  stall  power  are  calculated  and  subtracted 
from  the  residual  power  term.  At  the  difference  in  power 

(iiHP)  between  the  two  gross  weights  will  reflect  the  amount  of 
compressibility  power  through  the  AM  calculation  which  is  a 
function  of  tc ,( -M  = M>pip  - Mcr.+  .75  tc).  Now  a slope  of 
AHPAAM2  - AM^)  is  computed  where  1 and  2 refer  to  the  gross 
weights.  Then  values  of  M,-.r  are  swept  until  the  value  of 
(HP^o  - HPci)/(£M2  - AMj)  is  equal  to  the  computed  value  of 
AHP/t AM2  - AM^).  Knowing  the  value  of  Mcr,  the  compressibility 
power  is  calculated  and  subtracted  from  the  residual  power  term. 

The  value  of  tQ  is  determined  from  consideration  of  the  residual 
power  at  Vj^  power  speed  for  the  lighter  gross  weight. 

The  power  required  by  the  input  value  of  drajj  area  is  subtracted 
from  the  residual  power  at  the  speed  for  minimum  power  for  the 
lighter  gross  weight.  At  Vmax*  the  power  required  by  60  is 
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Figure  30.  Flow  Diagram  for  Input  Coefficient 
Program 
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subtracted  from  the  residual  power  term  leaving  the  power 
required  for  f.  The  equivalent  drag  is  calculated  from  this 
power.  Finally  the  total  power  required  is  calculated  based  on 
the  computed  coefficients  and  compared  to  the  input  value  of 
horsepower. 

Items  of  importance  to  note  when  using  the  above  method  are 
discussed  below.  The  value  of  the  weighting  factor  used  for  the 
induced  velocity  strongly  affects  the  value  of  £2-  the  values 
of  £2  are  unreasonable,  then  the  weighting  factor  K3  is  probably 
incorrect.  Some  typical  values  of  K3  are  given  in  Figure  31  for 
different  helicopters.  The  method  by  which  these  values  are 
chosen  follows. 

First,  the  power  is  computed  for  hover  using  standard  values  of 
input  coefficients  such  as  the  ones  for  the  AH-lG  helicopter. 

Then  the  value  of  K3  is  increased  until  correlation  is  improved. 
However,  reasonable  values  of  60 , 62,  and  Mcr  must  be  used. 

This  factor  is  used  to  increase  the  induced  velocity  to  improve 
power  correlation  in  the  low  speed  range  and  account  for  download 
on  the  fuselage  in  hover. 

Also  some  modification  is  required  in  the  angle  of  attack 
expression  for  the  different  helicopters.  This  change  yields 
improved  correlation  while  maintaining  reasonable  values  of 
£>o  and  £9.  The  average  angle  of  attack  of  the  rotor  is  approxi- 

/ 7gt\ 

mated  by  l-jjj-/.  However,  for  the  CH-3C  and  CH-53  helicopters, 
better  results  were  achieved  by  using  inPut 

coefficient  program,  the  afigle  of  attack  is  approximated  by 

Values  of  for  the  various  helicopters  are  presented 
in  Figure  31. 


The  only  parameter  not  directly  determined  by  the  INPUT  program 
is  the  amount  of  stall  power.  The  user  must  inpuc  the  divergent 
t~  versus  p curve.  The  value  of  f is  influenced  by  the  amount 
of  stall  power.  For  example,  if  a small  amount  of  stall  power 
is  calculated  at  Vj^^,  the  program  will  calculate  a large  value 
of  f to  match  the  power  at  that  point.  Ihis  f results  in  too 
much  power  in  the  immediate  speed  range  and  the  percent  error 
will  be  large.  By  increasing  the  amount  of  stall  power,  the  new 
f will  decrease  and  result  in  improved  correlation. 

To  aid  the  user  in  determining  the  proper  divergent  tc  versus  p 
curve,  data  from  Reference  18  are  presented  in  terms  of  a tc 
versus  p curve  in  Figure  32.  This  data  approximates  the  tc  value 
at  which  stall  is  encountered.  This  region  is  bounded  by  the 
lower  line  for  light  stall  (a  = 12  degrees)  and  the  upper  line 
(a  = 16  degrees)  for  heavy  stall.  The  divergent  tc  versus  p 
curve  from  Equation  47  is  also  shown. 
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HELICOPTER 

V 

V* 

AH-1G 

2.14 

7 

UH-1H 

2.30 

7 

CH-3C 

2.50 

9 

CH-53 

2.50 

9 

0H-6A 

2.14 

/ 

0H-58A 

2.14 

7 

★ 

★ * 


(t*-.14)  K3 


Figure  31.  Typical  Values  of  K3  and 

for  Various  Helicopters 
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Results 


Several  current  operational  helicopters  were  considered  in  the 
evaluation  of  the  equations.  The  necessary  physical  parameters 
for  these  helicopters  are  given  in  Table  I.  A set  of  input 
coefficients  foi  these  helicopters  is  presented  in  Table  II. 

A comparison  between  predicted  data  and  flight  test  data  for 
the  various  helicopters  is  presented  in  Figures  33  through  38. 

The  predicted  data  used  the  coefficients  presented  in  Table  II. 
This  comparison  was  made  using  representative  level  flight 
performance  for  the  different  helicopters.  Based  on  the  flight 
test  data  used,  the  equations  yield  satisfactory  results  for  use 
in  simulating  flight  paths  of  these  helicopters. 

Example  Problem 

The  AH-1G  helicopter  is  used  to  demonstrate  the  above  methodology. 
A sample  input/output  of  the  input  coefficient  program  (Appendix 
C)  is  presented  in  Figure  39.  Power  required  versus  airspeed 
is  usually  presented  as  shown  in  Figure  36  or  Figure  38.  If 
the  data  do  not  extend  to  include  hover,  then  hovering  out-of- 
ground-ef f ect  power  data  are  needed  as  is  the  case  in  Figure  36. 
The  data  of  Figure  36  can  be  reduced  to  dimensional  terms  using 
the  following  equations: 


GW 

(48) 

pA(SR)2 

HP  550 

(49) 

PA  UR) 3 

II 

(50) 

From  Figure  38.  the  power  required  data  for  the  AH-1G  helicopter 
for  two  different  gross  weights  can  be  obtained.  It  should  be 
noted  that  the  power  data  at  different  gross  weights  must  be  at 
the  same  flight  conditions,  i.e.,  rotor  rpm,  density,  and  Mach 
number . 

The  order  in  which  the  power  required  data  are  input  is  important. 
A one-to-one  correspondence  must  exist  between  the  power  required 
and  velocity.  In  addition,  the  first  point  must  always  be  the 
hover  point , while  the  second  point  is  the  power  required  at  the 
velocity  for  minimum  power  required.  The  input  value  of  m deter- 
mines which  input  velocity  (vc)  is  Vj^^.  In  the  example,  m = 7 
means  Vj^x  = 150  knots,  m = 6 would  mean  = 140  knots,  etc. 

Intermediate  velocities  are  used  to  give  insight  into  the  degree 
of  correlation  obtained. 


299-099-557 


55 


299-099-557 


56 


Dimensional  Data  for  Various  Helicopters 


Helicopter 

Conf  igurat  ion 

f 

'■1 

62 

M 

cr 

AH-1G 

Hog 

19.5 

.0075 

0 

1 

.75 

0H-6A 

Clean 

6 

.0075 

0 

.59 

.75 

UH-1H 

Clean 

20 

.0075 

0 

.59 

.77 

OH-58A 

Clean 

9.7 

.007  5 

0 

.59 

.75 

CH-3C 

Cl  ean 

35 

.0075 

0 

.59 

.75 

CH-53A 

Clean 

60 

.0085 

0 

.59 

.75 

Table  II.  Input  Coefficients  for  Various 
He licopt  ers 
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Figure  37.  Level  Flight  Performance  for 
HH-53C  Helicopter 
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Figure  39.  Example  of  Input  and  Output  from 
Input  Coefficient  Program 
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The  physical  dimensions  required  (Table  I)  are  rotor  radius 
(feet),  rotor  chord  (feet),  number  of  blades,  rotor  tip  speed 
( feet  per  second),  K3  and  K5 . Also  required  are  altitude  (feet) 
and  temperature  (degree  Centigrade)  for  the  given  data.  The 
lower  gross  weight  and  the  upper  gross  weight  are  the  final  data 
required . 

The  input  data  are  compared  with  the  c omputed  data  using  the  values 
of  f,  to.  ^1*  and  Mcr  determined  in  the  program.  The  actual 
error  and  percent  error  are  printed  for  each  speed  power  polar 
at  both  gross  weights.  The  values  of  f,  ho*  ^2*  and  Mcr 
determined  in  the  program  and  used  to  compute  horsepower  required 
are  printed  as  can  be  seen  in  Figure  39. 

The  sensitivity  of  each  of  the  parameters  as  related  to  magnitude 
can  be  seen  in  Figures  24  through  28.  Some  variation  from  the 
computed  values  may  be  used  with  little  change  in  the  power 
required.  This  is  the  reason  for  the  differences  in  the  values 
of  Table  II  for  the  AH-1G  and  the  value  calculated  by  the  INPUT 
coefficient  program.  Also  it  should  be  understood  that  a set 
of  coefficients  do  not  define  a unique  solution,  i .e . , there 
could  be  several  sets  of  coefficients  which  could  result  in 
approximately  the  same  degree  of  correlation. 
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APPLICATIONS 


Terrain  Following 

The  equations  developed  above  can  be  integrated  into  a terrain 
following  model  to  generate  an  accurate  flight  profile.  An 
example  of  such  a terrain  following  model  can  be  found  in 
Reference  25.  In  this  model,  the  aircraft's  terrain  following 
ability  can  be  investigated  as  a function  of  (1)  the  aircraft's 
characteristics,  (2)  terrain  following  subsystem,  (3)  load 
factor  limits,  (4)  command  clearance  altitude,  (5)  command  speed, 
(6)  terrain  type , and  other  operational  parameters.  The  pilot 
logic  was  developed  from  an  energy  viewpoint  along  with  some 
logical  decisions  about  how  the  pilot  flies  the  helicopter.  The 
pilot  logic  is  shown  in  Figure  40.  The  helicopter  will  always 
be  in  one  of  the  six  areas  of  speed  and  power  illustrated. in 
Figure  40.  Equation  1 (Figure  40)  was  derived  from  F = mV  and 
P = F-V  relationships.  It  is  used  to  determine  the  linear 
acceleration  (or  deceleration).  The  difference  between  the 
required  and  set  power  determine  (P.  Equation  2 (Figure  40) 
shows  the  use  of  the  weighting  function.  This  function  attempts 
to  recognize  that  pilots  do  not  correct  in  the  same  way  for  both 
large  and  small  deviations  from  the  command  speed.  The  weighting 
function  is  compared  to  unity  before  use  in  equation  2,  and 
the  smaller  value  is  used.  An  example  of  the  flight  profile 
generated  from  this  terrain  following  model  is  shown  in  Figure 
41. 

Decelerating  Turns 

The  question  of  penetration  distance  and  time  to  execute  a 180- 
degree  turn  at  constant  altitude  are  frequently  used  in  compari- 
son of  different  helicopters.  The  following  methodology  was 
developed  for  prediction  of  distance  and  time  for  a decelerating 
turn. 

The  pilot  can  choose  to  decelerate  at  constant  altitude  in  order 
to  supply  additional  power  for  a high-g  turn.  The  time  history 
of  the  flight  path  must  be  constructed  of  a series  of  circular 
arcs  since  speed,  radius,  and  g- level  are  varying  throughout  the 
turn.  The  procedure  is  described  below. 

(1)  Determine  the  variation  with  speed  of  the  maximum 
transient  thrust  capability  of  the  rotor  (see 
Maneuverability  Limits).  Divide  this  by  the  weight 
to  get  the  V-n  curve.  This  maximum  thrust  capability 
may  be  limited  by  rotor  instabilities  like  flutter  or 
weaving,  or  by  loads,  vibration,  or  stall  --  depending 
on  the  specific  rotor  system. 

(2)  Determine  the  power  required  at  the  maximum  rotor 
thrust  using  the  power  equations. 
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(3)  Determine  how  deceleration  varies  with  speed  using 
equation  39. 

(4)  Select  the  time  interval  to  be  used.  One  second  is 
usually  sufficient. 

(5)  Determine,  the  average  deceleration  for  the  first  time 
interval.  Then  compute  the  average  speed  and  determine 
the  radius  of  turn  from 

R = V2/fg  '/n2  - l) 

(6)  The  amount  of  arc  traversed  in  the  first  time  interval 
is  given  by  d./dt  = V/R  for  a ft  of  one  second.  If 
computed  instead  of  graphically  constructed,  the 
penetration  distance  is  calculated  as  follows: 


penetration  distance 


V cos 


heading  change' 


(7)  For  the  next  time  point  the  entry  speed  is  the  exit 
speed  of  the  first  time  interval.  Return  to  (5)  and 
continue  until  the  desired  heading  change  is  achieved. 

An  example  of  a 180-degrce  decelerating  turn  is  shown  in  Figure 
42.  A comparison  of  predicted  data  using  the  above  procedure 
and  measured  data  from  Reference  7 for  the  AH-1G  helicopter  is 
presented  in  Figure  42.  Entry  airspeed  for  the  above  example 
was  150  knots. 


Low  Speed  Maneuvering 


The  helicopter  has  a capability  of  accelerating  in  any  direction 
from  hovering  flight.  This  can  be  used  to  increase  the  surviva- 
bility at  very  low  speeds  since  changing  the  direction  and 
magnitude  of  flight  path  acceleration  complicates  the  ground 
gunner's  probLem.  Using  the  methodology  from  preceding  sections 
and  observing  physical  laws,  a method  (Reference  26)  for  deter- 
mining the  capability  of  the  helicopter  to  accelerate  both 
vertically  and  horizontally  at  the  same  time  is  presented  . In  high 
speed  flight  at  moderate  altitudes,  maximum  g-capab ility  can  be 
used  since  the  power  required  can  be  supplied  by  loss  of  energy 
(altitude  or  airspeed).  However,  rotor  maximum  g-capability 
cannot  be  used  effectively  in  low  speed  flight  since  kinetic, 
potential  and  rotor  energy  levels  are  low  and  nearly  constant. 
Therefore,  low  speed  maneuvers  are  done  at  constant  power. 
Horizontal  acceleration  is  most  efficiently  done  by  tilting  the 
rotor  thrust  vector  so  that  the  horizontal  component  acts  to 
accelerate  the  helicopter  as  illustrated  in  Figure  43.  The  main 
rotor  thrust  must  be  increased  to  maintain  altitude.  This 
requires  an  excess  of  power  above  that  required  to  hover  and  is 
the  means  by  which  excess  power  is  converted  to  accelerative 
force . 
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Figure  U2  - 180  Degree  Turning  Performance  far  the 
AH-1G  Helicopter 
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Figure  43.  Thrust  Vectoring  for  Linear  Acceleration. 
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The  power  required  by  a typical  high-performance  helicopter 
in  forward  flight  is  illustrated  in  Figure  44.  The  effects  of 
normal  acceleration  on  power  are  indicated  by  the  lines  of  con- 
stant g.  The  pilot  may  elect  to  use  the  normal  acceleration 
capability  to  accelerate  vertically  (increase  rate  of  climb  or 
arrest  rate  of  descent)  or  to  accelerate  parallel  to  the  ground 
while  maintaining  a 1-g  vertical  component  of  the  thrust  vector. 
The  power-limited  g-capability  may  be  determined  from  the  data  of 
Figure  44  by  noting  the  g-capabi 1 it y at  each  speed  at  a constant 
power  level.  Turning  flight  requires  acceleration  normal  to 
the  flight  path  and  the  power  required  that  is  associated  with 
producing  increased  thrust.  Level  flight  acce lerat ion _ requires 
power  (equation  39).  Thus,  the  maximum  acceleration,  V,  is 
determined  by  the  excess  power  available  and  the  velocity.  It 
should  be  noted  that  the  predicted^  tends  to  infinity  as  V tends 
to  0 which  is  a limitation  of  the  V equation.  A graphical 
presentation  of  the  V equation  is  given  in  Figure  45  for  a 
helicopter  of  12,500  pound.1  gross  weight. 

The  preceding  discussion  has  assumed  that  there  is  no  gain  or 
loss  of  altitude  (potential  energy).  The  power  required  to 
climb  is  given  by  equation  28.  To  account  for  climb  and  descent 
effects,  the  power  computed  by  equation  28  is  .-.dded  to  (descent) 
or  subtracted  from  (climb)  the  maximum  power  available.  A 
horizontal  line  is  drawn  across  Figure  44  at  the  new  power  avail- 
able and  the  corresponding  acceleration  capabilities  are 
determined  as  before.  For  example,  a 10  ft/sec  rate  of  climb 
would  require  267.4  horsepower  for  the  example  helicopter.  The 
maximum  power  available  then  would  be  1830  - 267.4  = 1562.6 
horsepower.  The  normal  acceleration  in  hover  would  then  decrease 
from  1.18  to  about  1.05  in  a 600  ft/min  vertical  climb. 

The  power-limited  linear  (flight  path)  acceleration  at  a specified 
value  of  normal  acceleration  can  be  determined  from  the  data  of 
Figures  44  and  45.  For  example,  if  the  helicopter  is  banked  45°, 
sustaining  1.414  g in  level  flight,  the  power  required  at  60  knots 
is  1110  horsepower  from  Figure  44.  This  gives  an  excess  of 
720  horsepower  which  may  be  used  to  increase  flight  path  velocity. 
Entering  F:gure  45  with  720  horsepower  at  60  knots  gives  about 
0.31  g (5.5  kt/sec)  linear  acceleration.  This  normal  and  linear 
acceleration  is  approximately  correct  except  that  the  vector  sum 
of  1.414  and  0.31  is  1.447.  Since  the  power  required  data  of 
Figure  44  refers  to  the  total  thrust  produced,  the  vector  sum  of 
the  normal  and  linear  acceleration  must  be  used  to  obtain  the 
power  required.  This  involves  reentering  Figure  44  at  the  higher 
normal  acceleration  (1.447)  and  obtaining  the  power  required. 
Figure  45  then  shows  a linear  acceleration  attainable  of  0.29. 

This  iterative  process  continues  until  accelerations  are  deter- 
mined which  are  compatible  with  power  and  vector  summation 
constraints.  Figure  46  is  a combination  of  Figures  44  and  45 
which  indicates  the  approximate  combined  acceleration  capability 
if  the  vector  sum  is  not  included. 
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Figure  45.  Linear  Acceleration  Produced  by  Excess  Power 


Figure  46.  Approximate  Combined  Normal  and  Linear  Acceleration  Capability 
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To  determine  accurately  the  combined  acceleration  capability, 
the  data  of  Figure  46  are  replotted  in  Figure  47  at  speeds  of 
interest.  The  curves  starting  at  the  1.0,  1.2,  1.4  and  1.6  g 
values  of  n are  constructed  from  the  geometric  relationships  of 
a right  triangle  as  indicated.  The  values  of  n()  and  n^  at 
the  points  where  the  parabolas  intersect  the  speed  lines  obtained 
from  Figure  46  satisfy  both  power  and  geometric  constraints. 

For  example,  from  Figure  46  the  indicated  linear  acceleration 
capability  in  level  flight  (n=1.0)  at  30  knots  is  0.74g.  This 
is  point  A on  Figure  47.  The  rest  of  the  30-knot  line  of  Figure 
47  is  determined  by  reading  corresponding  values  of  n and  nl( 
(linear  acceleration)  from  Figure  46  along  the  vertical  30-knot 
line.  Where  the  30-knot  line  intersects  the  geometric  curves 
previously  described,  both  power  and  geometric  constraints  are 
satisfied.  Point  B on  Figure  47  indicates  a level  flight 
acceleration  of  0.56  g at  30  knots.  The  normal  and  linear 
acceleration  capability  for  the  example  helicopter  is  given  in 
Figure  48. 


In  hovering  flight  the  above  method  fails  since  the  linear 
acceleration  equation  is  infinite.  The  maneuvering  limits  near 
hover  are  determined  by  the  horizontal  and  vertical  acceleration 
which  results  from  directing  a thrust  vector  of  constant  length, 
1.18  g for  the  example  case.  The  magnitude  of  this  thrust  is 
determined  by  power  available  as  indicated  in  Figure  44.  The 
result  is  a circular  arc  with  a radius  equal  to  the  g capability 
with  V = 0.  The  above  combined  acceleration  capability  is 
graphically  portrayed  in  Figure  49  for  low  speed  flight.  Figure 
50  is  a similar  presentation  which  shows  data  from  hover  to 
maximum  speed.  The  kinetic  energy  available  at  high  speeds 
enables  considerably  higher  g-levcls  to  be  sustained  for  short 
periods  to  make  rapid  changes  in  the  flight  path.  The  power 
required  for  these  maneuvers  is  provided  by  loss  of  speed 
(kinetic  energy)  during  the  maneuver.  Altitude  (potential 
energy)  can  be  lost  to  gain  power  tor  maneuvering  at  g-levels 
outside  the  envelope  indicated  in  Figure  50. 

At  this  time,  no  flight  test  data  are  available  to  validate  this 
theory . 

Survivability 

Survivability  studies  require  some  type  of  flight  path  deter- 
mination. The  impact  of  the  accuracy  of  these  trajectories 
varies  depending  on  the  survivability  model  itself.  This 
simplified  method  for  determination  of  flight  trajectories  offers 
accurate  representation  of  the  flight  path  trajectories  for 
minimum  computer  run  time  and  computer  storage  space. 
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Figure  49.  Effect  of  Speed  on  Normal  and  Linear  Acceleration  at  Constant  Power. 
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Flight  Test 

The  evaluation  of  the  maneuver  capability  of  a helicopter  in 
flight  test  often  results  in  data  which  reflect  pilot  technique 
more  than  the  actual  capability  of  the  helicopter.  From 
Reference  4,  return-to-target  maneuvers  were  performed  by  three 
different  pilots  using  the  AH-1G  helicopter.  A standard  devia- 
tion of  1.01  seconds  was  calculated  for  level  turns.  From  a 
statistical  viewpoint,  this  is  a very  large  scatter  band.  The 
large  variations  observed  in  the  test  results  degrade  precise 
quantitative  comparisons  of  maneuvering  performance.  If  one  is 
faced  with  comparing  the  maneuver  capability  of  two  different 
helicopters,  the  influence  of  pilot  technique  may  reduce  or 
augment  the  helicopter’s  performance.  For  this  reason,  a method 
is  needed  for  the  evaluator  to  compare  maneuver  capability  of 
different  helicopters  independent  of  pilot  technique.  The 
analytical  methods  to  predict  acceleration,  deceleration,  and 
turning  performance  are  presented  in  the  preceding  sections. 

These  parameters,  combined  with  conventional  performance 
parameters,  define  the  total  performance  capability  of  the 
helicopter.  Almost  all  maneuvers  are  cor  binations  of  one  or  more 
of  the  above  basic  maneuvers,  and  energy  methods  may  be  used  to 
combine  these  basic  maneuvers.  The  total  performance  capability 
of  different  helicopters  may  be  quickly  compared  by  the  overlay 
of  the  energy  diagram  such  as  Figures  3 and  4.  In  addition, 
since  power  required  at  different  altitudes  is  calculated  by  the 
power  equations  for  the  energy  diagram,  the  fuel  flow  can  be 
determined  depending  on  the  engine.  Then  a chart  of  specific 
range  and  endurance  may  be  constructed  as  a function  of  altitude 
and  airspeed  as  shown  in  Reference  27.  To  nroperly  consider 
differences  in  gross  weights  and  loading  conditions  with  differ- 
ent helicopters,  it  would  appear  reasonable  to  compare  total 
performance  at  a gross  weight  where  both  configurations  can 
exert  the  same  maximum  ”g"  load.  For  example,  consider  a compari- 
son of  the  charts  at  a ''g’’  load  of  1.5.  The  energy  diagram  for 
the  AH-1G  helicopter  at  a ’’g"  load  of  1.5  is  given  in  Figure  5. 
This  would  assure  the  evaluator  ot  a given  level  of  maneuvera- 
bility while  comparing  other  aspects  of  performance.  Since  the 
techniques  used  to  develop  g-capabili tv  are  primarily  based  on 
energy  and  momentum  considerations,  the  efforts  of  different 
configurations  on  performance  parameters  are.  minimized. 
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RECOM*  ’ENDAT 1 ON’S 

1.  Additional  work  should  be  initiated  lor  continued  develop- 
ment ot  the  low  speed  maneuvering  theory. 

2.  The  existing  theory  should  be  validated  with  flight  test 
data  and  the  rotor  degree  of  1 reedoin  added  to  the  equations 
tor  low  s pe e d . 

3.  Methods  ot  improving  the  helicopter's  capability  in  the  low- 
speed  regime  should  be  investigated. 

4.  Investigate  further  the  effects  of  blade  flapping  inertia  on 
the  maximum  tc  versus  relationship. 
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DERIVATION  OF  BANK  ANGLE.  LDAP  FACTOR  RELATIONSHIP 

Equation  (7)  was  derived  from  equations  on  page  116  of  Reference 
17  as  follows: 


Assume. the  aircraft  is  turning  about  a vertical  axis.  Thus, 
Euler  and  ; are  zero  and  the  body  angular  rates  become 


P 

q 

r 


where 

p = roll  velocity 
q = pitch  velocity 
r = yaw  ve loc i ty 


- . s in- 

v cos  » sin;  (A-l) 

• cos  e o s •+ 


= pitch  attitude 
; = roll  attitude 
t - yaw  attitude 


These  angular  rates  are  then  substituted  into  the  generalized 
equations  of  force  equi 1 ibri um  : 


X 

Y 

Z 


- mg  sin 
+ mg  cos-. 

+ mg  cos 


= ml  u + qw  - ru) 
sin:  = m ( v ♦ ru  - pw ) 
cos  ; = m(w  + pv  - qu) 


(A-2 ) 


where 


X 

Y 

Z 

m 

u 

v 

w 


force  on 
force  on 
force  on 
a i rc  ra  f t 
ve locity 
velocity 
ve locity 


vehicle  in  X direction 
vehicle  in  Y direction 
vehicle  in  Z direction 
mass 

in  X direction 
in  Y direction 
in  Z direction 


Assume,  for  coordinated  flight,  that.Y  = 0.  Also  assume  no 
change  in  linear  velocities  so  that  u,  v,  w are  zero.  Thus, 

g cos&  sinf  = ru  - pw  = u-(cos“  cosC + sin&  tana)  (A-3) 


if  we  define 


tana 


w 

u 


and 


tanp  = J 
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Now  solving  equation  A-3  for  — 

Q 


jjv  _ tan ; 
g " 1 ♦ K 


^ _ tani-  tana 


(A-4) 


Use  equation  A-l  in  the  Z force  equation  and  substitute  eqiation 
A-4  for  u*/g.  Equate  Z/mg  to  -n: 


n = ,siH’ — tarw  _ u ...  cos_sin2  . cogf,  cos4,  (A.5) 

g ? 


Then , 


n cos.  = cos 


2.  . , sin  . . sin;  tani  sin(-  , 

cos  . ♦ cos  * t — — rr  i — r— r CA-t 


Consider  the  following: 

1.  If  a = (»  = *»  = 0,  then  the  usual  equation  is  obtained. 

n cos ; = 1 

2.  If  a = t - 0,  then  1 + K is  unity  and  the  third  term  of 
equation  A-6  is  zero. 


Since 


2 2 
sin  : + cos  t = 1 


n cos>  = cos- 
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ENERGY  DIAGRAM  PROGRAM 

/•  ENERGY  KANELWEHAUILl  TY  HV  DIAGRAM  • /; 

/•  Innut  stcilun*/; 

_ ?UT  UiTI'lf  .1>J1  <nd  P>^.  del  lln«  H i(w).in(n)  ')■ 

CEf  llSI(.i.n); 

k*n*l ; 
w*n*  1 ; 

UECLAMt  a ( 21 ) , xnUU)  ; 

I f .tl  lucjtf ) THEN  FHEE  f,ran; 

_AUOCATf  ntu^rMidOl 

DECLARE  f DEC ( b ) CONf ROIL ED  , rin  DECI I ) CONTROLLED  ; 

DECLARE  a»5U1u  EXT  ENTRY; 

GET  Ll  jTIrt^^chorl.bladt,.  t»n<Lrlrae.  w lah. If ); 

GET  L I :.T(uat,eta) 

GET  LI‘.»T(vlnc,vlnl  t , h I nc, h I n 1 1 ) ; 

(•  .ir^'cijUtlon  n»  <t  util  »r Id  »olnt«/j 

Bo  j-I  To  m*1; 

IJU  1-1  to  ml; 

h«<  1-1  )•■•! fief'll n 1 1;  _ ____________ 

■l»td*(  1- , uuuujttl  Vj  35*h)**5. 2Sb"l; 
tsl  ao4t*2M,lb; 

- »!rl#uil"^wmz^l2giu 

• I &na'l«  td/ theta; 
t • theta* t*  1 ; 

. • ii tCiJi  ... . 

drha*s  l cii*.  JU2  J ; 
ve1*vlnc*j*«lnl Inc; 

.(All  d e 3 wlf UA^-UmiiULUAA^ Ufld»rirAA»anlAh,xAl,lf,rlrhn,  rina,dh»h 

IF  T»lf M GO  TO  c5U; 

tWM'lluu; 

GO  Tj  cU;  _ 

h.m  — .<i2*(  • -VJUU)*1100; 

,»*•  >*>!)•(  .pa-dhp)/  (eta*welgh) ; 

. 11  * aJi  VU 

EMU  «.i; 

I Mn  c2; 

/•  find  m » I nun  valui  of  oi  AL  tich  jralnell*  I nrraaiant  t/ j 

30  j*l  TO  n*  1 ; 
iu  1-1  TO  »»♦!; 

_*dJ-f(ULLi 

END  c/; 

CALL  order(A,'i); 

ran{ 1, j ) *a(m);  _ _ 

ra. i(2,j ) • a ( 1 ) ; 

END  Cl; 

_ L • ..deter  I n|t  |yn  Bf  LlflJUL  Of  OXLU 

00  j«l  To  n; 

«n( j )«ran( l,j  J ; 

CNU  clU;  . „ ... 

CALL  order ( *n,n J; 

*i*i  >ax*«n(  n)  ; 

- aha* ■; rune UiijAAi^AUiiLaiaajLlU  — , 

DO  j«l  TO  »>; 

*n( j } *ran(2, J ); 

EMU  Cl  TJi  _ ........ 

CALL  order ( *n,n) ; 
pt<iln>jtn(  1); 

_ A'J  nnrunclumlnt.S+iUnt  plain)  1; ■ — 

nnax»l  runc( Ana*/ i )*i ; 
mil  n*  t rune  ( xml n/ 5 )*5 ; 

#n ■{ xnax-xnl n)/5;  - 

an* an* 1; 

PUT  L I ST ( an, xmax, xnl n) ; 

IF  aHoca(psc)  THEN  FREE  p»c#vc,hc; 

ALLOCATE  ps c ( an ), vc( an, n ) , hc( an,n) ; 

PECLAKE  n* c DEC(6)  CONTROLLED  , vc  DEC ( 6 } CONTROLLED  , he  DEC(I)  CONTROLLED  ; 

psc-Q;  - 

vc,hc-0; 
p*c(  1)-a<  tax; 

. PO  lb- 2 TO  an;  — 

p»c(  i b ) *p*c  ( I6-1)-S; 

|NU  Clb; 

/•  p*c  sweep*/;  ._  .... 

DO  lb*l  TO  an; 

/•  velocity  grid  sweep*/; 

_i>0  j-1  TO  

1*1; 

IF  rand,  J ) >*p»c(  16)  THEN  GO  TO  dl; 

IF  ran(l,j*n)*U»c(l6)  THEM  GO  IQill*  - - 

60  TO  cl  8; 

IF  rant  l,j  *l)>psc(  16)  THEN  GO  TO  <13; 

vc(  lb,J  )«vlnc*(J*U*ylailr\£lacj 

lie ( I 6,J  ) *0; 
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CO  TO  dl; 

yc(  I b,J  )■»  I ncMuictl  fcLf  ( » - 

lic(  lb#}  )*0; 

CO  TO  ell; 

_|F  ^nUJj^dUJ  1HU  6Q  Ifl  dki 

*c<  I 6, j )»vl nc*J *vl nl t-vlnc; 
hc( I b, j )*0; 

60  TO  ell; 

00  1*2  TO  m; 

I F l,J  X-p»c(l6)  THEM  GO  TO  dS; 

. JJO  To  

If  f<  I,  i Xptc(lb)  THEN  GO  TO  df; 

Y c( lb#])>vlnc*j«vlnlt*vlnc; 

00  To  ell; 

YC(  I to,  J )>vlnc«j*vln|  fvlnc; 

- iic  c i 

00  To  cl»; 

|N0  iU; 

(NO  al ; „ 

(Nl)  cl/; 

00  U-i  To  fcn; 

~tto  j-l  jy^-L 

|f  osc  ( lb)<«f(l,j)Aptc(  l6)>»Ml,|*i)  THEN  CO  TO  cSO; 

00  TO  cSl; 

yc(  lb#JM)«v|nc«(fll#Jl-flAc(lbn/lfll#J)-Ul,JtUiaaUcala*i*44«*m#A 

00  TO  c j 4 ; 

Hie  ' lfc)>«Ml,j  M|,«C(  I6)<*f(  J,J  »1)  Then  00  TO  c52j 

_Hu  Tg  

»e(  I6.J  )•*!«•(  p»e(  It ) -Ml,  J ) )/  ( t(  1,  J • 1 )-f  ( 1,  J > )»»lnc*J*»lnU-»lncj 

(NO  r.  i»  b ; 

(NO  ciG; 

/•  output*/. 

PUT  LUT(  ' * ); 

_ PO  T Will' a*  * >j — 

flu  1 1<  • X To  am; 

DO  j»l  Tj  n; 

|F  vc(  U,j )»0‘«htULil*O  IHIM  6Q  TO  Ul;  . _ _ 

*UT  IMAGE  (.>SC<  lt),ve(l6,J>,hc(li,J))<ln2); 

IMAGE; 


ell:  (MU  c2  5; 

(NO  c2to; 

/•  * »ubrouLlaA,Arr^a«AA  a caluan.ArrA*  Iaao  «aIm«  •#  itttndlne  ir<ir»Ji 

ordtri  PHOCEUUHE  (d*,dm);  ■■■■"■in*  wrmwr^f^ 

bit  00  12*1  TO  dn; 

I S • I 2* 1; 

bJ:  00  dl -IS  TO  di; 

IF  d*  ( 1 2 )<»d*(d! ) THEN  GO  TO  b4; 

(— - — 

4*( I 2)«d«(d1 ); 
d«(dl )"t«np; 

_nki (HQ  L\; 

(HO  66; 

(NO  order; 

•t)010l  PMUCCOUHE  < r, c, b,wr, f I , gwl , vl,  nl , rho, vi,hp) ; 

/•  THIS  SUBROUTINE  CALCULATES  POWER  REQUIRED  FOR  ENERGY  MANEUVERABILITY  •/; 

f constants  for  wprwl  HeHcoptcts  il  m ttft  *H 

list  de|U*.0075; 

de  j 7»1 ; 

• •(.21;  - - - - 

ncro*. 7S; 

k 1 - 1 ; 



/•  compute  tome  constants  •/; 

SO:  p I ■). 1 4 1 5 9265  55; 

«raa-f»l*r«*2j  — 

»l|*b*c/(plar); 
nondf • rho*area*wr*«2; 


CtQ*7/(slg*a); 

f-fl; 

«w*tfwlj  

v-y1  *1.  f>8  70; 

mu»v/wr ; 

, HHt»} 

IF  nu>. 14  THEM  kl-1; 

ny2«i  hj**2 ; 

*3ewr«(  J*mu)/vs;  

n«nl ; 

/•  compute  some  more  constants  */; 

; 

ct •n*gw/nondf ; 

Ctn»snrt(*4*x4/ (nondf •nondf )*ct#ct); 

« 1 •*4«v/ ( ctp*nondf ) ; 


T 
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«2*ctP*noiw1f  / ( 2*rHoaar«a) ; 

/*  compute  Induced  velocity  from  momentum  theory  equations  */j 

eownf  #; 

I F v>i»2.2  THEN  CO  TO  »17; 
vl  *1. 0W*»qrtU2) ; 

IF  v-0  THEN  » TO  ill;  

• 17i  vln*»2/v; 

ill  vl«al*«2/  <v*v*,IM*vln*vln)**,$j 

+»  vImH.MI  THEM  00  TO  i!l; 

count ■count* I j 

IF  count  > kQ  THEN  CO  TO  '.O; 

vln-,S*v»  ♦ . k*vlnj  ... — 

GO  TU  *9; 

ilOt  FUT  L I ST( * ca  1 culat Ion  non-convergent, I nduced  velocity  to  0*)j 

- 

GO  Tu  *12; 

• Ut  vl«vl-*l; 

/•  compute  oowtr  required  tvcTton  Vr  

/*  flat  plata  drag  power  */; 

•12 i h0l*xt*v/5SU; 

1*  ■ Indore*  now? — *77  ■ ■ ■ ~ 

hp2*kl*ctp*nondf*vl  / S SU; 

/•  bled*  profile  drag  power  at  taro  lift  •/ j 

np|*delu*alc*nonftp*Tl*t.i*matTAtT — — 

/•  blade  profile  drag  power  caused  by  ancle  of  attack  •/; 
tc *2 *ctp/ * I s ; 

trcrr*«lw.t/»nf  tflvOOmun  j 

hp**0; 

hpk-del 2*alf *nondp«( l*fc.fi*mu2)*(ct*ctc)**2/l; 

IF  tc< tccri  THEM  00  TO  111;  * - ‘ “ 

del tc*tc-tccrt ; 
hp|*200u0*deltc**l,S; 

/•  compressibility  power  •/; 

• Hi  IF  «)>fKro  THEN  GO  TO  si); 

dcp-U; 

00  TO  Ilk  J " - ’ 

• Hi  de1it*a)*ncro*  ,7S«(2*Ctp/sl|); 

dco*delm**S*<  .00JS-dslm*(.tt22-.ll*dalm)); 

alb  i l'ipO»kf*ifvu*Hmwt|»; 

bp«hpl*hp2*hp3*hpS*hp$*hp» ; 

END  aaJOlO; 
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MANEUVER  DIAGRAM  PROGRAM 


r 2 i 

cl  t 


clO: 


c 9 : 


clO: 


Cl  5 : 


cl6 1 


clU 


a 1: 


/•  MAH!;.  RAblllTV  MANf  IIVF  H MAGMA*  •/; 

4 • 1 n*Jul  atutluu*/.  - _ 

P T t I .T  ( 1 I ♦ *»  il  and  i »i»U,  dr  I I I n#  19  a (m), *n( n ) '); 

git  n ,n 

k-  i*  1 ; 

w»*l; 

.-I  Cl  A-.  l «(  :i  > , »-.Uu)  ; 

li  alluwUll  liti-it  liii-1  J +AdAi — 

ALL  (‘ATI  * * fc,  *),  ranl^Nl ; 

JtClA  I ♦ Kit*)  CUNT  do  l LEU  , ran  UECU)  CONTROLLED  ; 
uEClAUI  as  lulU  UT  ENTHT. 

GET  L 1 j T < rad,  c*»ur  t,  r>1ad«,  t md,  drag , wt  I gh,  h) ; 
f,C  T l I . T i .<<» t,  eta  ) ; 

. £i  T LUU  < li«U.»aUtC»aiJ>fcU  ... 

/•  r.ii  r;  , ' ,i{  I o>’  of  ,»s  at  aacN  grid  uolnt«/j 

l • » T ; >*l; 

■>w  I'l  Tj  >*«)  — . 

If I • l )•  .1  ,.r; 

1st  •( 1-..  I...  .t»T  J .•n#»«5.25U; 

- - Aalaw«t»^«'A-*W  ...» 

t'  *•» 

s • j • • : s t * / l ••  t a ; 

l»l  » 

! * > ■ *j  •• , n 1 ' Si.ij  *t'ir  t ( t ) , 

Ir  ■ i ■ * s 1 1,  • • ■ U v * W A : 

»«i*<  Inci  > — _ . , . — . — 

C A l l .is  l , 1 , ( r.»  1,  c'mrd,  h 1 ad#,  tmd,  dr  ag,  w*  I gh,  v*1 , I f , dr  Go, dvi,  dho) ; 

IE  T ' i f n GO  TO  cJO; 

■ >oaallw«,  - — - — — - — — 

’iJ  T K ■ 

■-..a*-.  *.’•  - 1 /00)*1100; 

dJipi;  I4UI  j 

Mi,,)#  >%. 

in 

INw  c*; 

/•  * I • * *•!  ur  ralue  of  pt  at  each  velocity  lncrtn#nt*/j 

...  •:  r , . 1. 

— la.  la  «a*l*  — 

*<  I ) - M I , j ) ; 
r *1  : 1 , 

CALL  urdarta,.-*!; 
rar  ( 1 , ; )•«(  ); 

r.»f’  (?,.)•»(  1 ’ ; 

- 4**w  *.ii  - • - — — ■ - 

/•  Jeter  iriatlo*'  of  rang#  of  ns#/; 

) ,•!  Tn  . 

«•«(  j ) >ran(  l,  j # » 

IN.,  rjy; 

CALI  if  ler  ( ; 

******  a <w>  4 «♦>»  - — ■ — ■ 

«•  <a*  #t  r.,Mc  ( •>%•'**♦  . '»#s  I gn(  nsna* ) ) ; 

no  ; ■ 1 TO  n; 

«n(  I ) »rar»(  2,  j ) i 
EV-  cl  n ; 

CALL  or da r(»n#n); 

- f*V*-*i  • ■ 

4'  I n«t  r jnr.  ( ;>»'■> I n*  , I gn{  i»sr  I n ) ) ; 

#'ia*«trMnr  ( * »a */ 0 ) •!> ; 

*n|n»trur>cl  *r>l  n/  S > • k ) 

an • ( *' ..1  *•*  1 1 n)/5; 

1 ; 

PUT  L I *»T(  .in,  n'  >a*,  anln ) ; 

IF  pic ) THEN  FHft  P»C,vc,nc; 

ALLOCATE  i>sc(dn),wc(an,n),nc<an,n); 

DECLARE  n»r  UfC<4>  COWTROlttO  , v«  UEC<A*  MHJHQ iiCO  , COHIftOUXO  * 

osc»U; 

VC, nc«U; 

- 1 )»w>«Ma» — — - 

00  lii "2  TO  an; 
p*c(  I ♦»)•.>*  C C I &-1  )“ij 
€NI>  cll»J 
/•  psc  sweep*/; 

DU  I fj - 1 TO  an; 

- / * velfseltr  — ■■  ■-  — - 

00  j -1  TO  n; 

1-1; 

IF  rar(  1,  j >>-rnic(  lfc>  THEN  GO  TO  HU 
|F  rant  1 , j ♦ 1 > > -PSC<  I b>  THEN  GO  TO  ri  2 J 

GO  T',  cl  a; 

♦ F rnnil,  THEN  GO  TO  4A» 

vc( IU,j ) • v I nr  * ( j ♦ 1 ) ♦ w i n I t-vlnc; 
nc  ( 1 1>,  | ) ■ 1 ; 


I 


299-099-557 


91 


BELL  HELICOPTER  commnv 


dll 


dll 


dbi 


dll 


H*T“ 

ill 
cl  1 1 

-rt*r 

CS5 ! 

dill 

dli 

nti 

ci*i 


cl*  i 
dll 


lull 


M 10  Cldf  - 

veil*,  I )**lnc«(o»c(  ll>*f(l,JI)/tr<l*l,J»l)*M  l,J>>»*lnc»J»»lBU-»lnc» 
net ld,J )«1| 



IF  rjnll,  J )>pdd  111  THIN  00  TO  d*i 
vcl  Id, J Idvlnedjdvlnl t-*lnci 

ncl  l*,J  »*1>  - • - 

00  TO  cllj 
00  1*1  TO  nj 

-IN-WfJ  I TWO  M TO  dli  

00  TO  «1) 

IF  f 1 1 , j Hotel  II)  THEN  00  TO  d*i 

»e(  I*,  j 

ncl  l*,])«l*<  I -1  Xnlnei 
GO  TO  eld) 

-rcl  l*/l  >»*ln«»l»oln**'^>)**)  1 1 1 

ne!  lt>,J  >*'it»e«cf  1 l-l,J>*Mc(H)>/(TCI-l,J)-Tll,J  ))»nlne»(  1*1 1*1) 

00  TO  dll 

(NO  d*j  - • 

END  *1| 

END  *'Tj 

m i*«i  to  nr  — — 

no  J*1  TO  n-ll 

IF  i»c(lt)<-f(l,j)IIPtc(U)»r(Lj*ll  THEN  OO  TO  cSO; 

OH  Tn  dli 

»e<  IH, J d ) •*ln««( f ( 1, J l-iiel  H))/(F(1,  J )-F(l,J»l) )•» lne»J»»lnl  t-vlnej 
iit(ll..J*ll*l) 


IF  Iitcl  I b I >*f 1 1, j I Node  ( HX*f(l,J*l)  THEN  60  T0  e»I» 

CO  Tf)  cib; 

»c(  16, J l-rtncdlnpcl  l»-P<Lf  1— pJppUUpnImi- 

ncl  I J Xlj 
ENU  ell; 

■tno  

/•  outiiut*/; 

PUT  LIST l*  lid  «c  ne  thd  pdthd  ’ll 


OO  l««l  TO  an;  - 

OO  j«l  TO  n; 

IF  ,cl  Ib.J  l-OAncl  14, J ><•!  THEN  SO  TO  Clll 

WllilildUlmlllill’MlHijVUdlMlldll  Ui))d|i>IH>| 

pithri-nicl  H>)»tlid/lJ 7. 5dll.ll) 

IF  ncl lb, J )•!  THEN  IhddOi 

PUT  image  1 peel  J I#  tN4»ppdMl<  Iwllp 

IMAGE ; 


Clll  EHO  Clli 

ENU  c Id; 

GET  llbTlt.tlnc.tlnlt); 

III  00  '-1  TO  Ml;  

«lid"t  I nc«(  I -1  X t Ini  t ; 
b/i  DO  j*l  TO  n.l; 

PIT  «vtHe»T)»TT»PTntti 

lfd|(tMdveld.000*15)ddldl>dd.S) 
PUT  IHACCIthd.vpI.IMIlnlll 
T«l|  IMAGE: 


• •• • • • • * 

t hu  hi; 


blO| 

b9t 

mu  ni; 

uu  i-i  to  i*«l; 

00  J-l  to  n*l; 

i f • !*( I -1 )*ntnc j 
vel-vlnc-(J*l)-vlnlt; 
thd-1092. l-l If ; 

"Ml  TWUCrfTf, ti'lbli 
IMAGE; 

CUD  b9; 

f NO  blU;  . . # 

/•  this  subroutine  arranges  e column  errey  Into  veluee  of 

orn#r  s 
b$ 1 

DO  12-1  TO  dm; 
15-12*11 

b J j 

TJO  JMT  lu  ta; 

IF  dilUX'dAldl)  THEN  00  TO  Ml 
ttmp-d«l 111; 

bit 

Oil  1 1)  )I 

d«(dl )*t*np ; 
ENO  bl; 

END  order; 
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SIMPLIFIED  THRUST  MODEL 


Below  is  a derivation  of  the  equations  used  in  the  simplified 
math  model. 


The  thrust  for  a given  number  of  blades  (b)  can  be  calculated  by 

_R  .R 

L O 

(B-l ) 


T = J bqCL  cdr  J | ^n2CL  cdr 


where  Cl  is  the  section  lift  coefficient  and  Vn  is  the  velocity 
normal  to  that  section.  The  velocity  normal  to  the  blade  is 


*Rx  ♦ V sin  • = (x  ♦ ,1  sin  • ) R 


(B-2) 


where  x = r R and  r is  the  distance  out  the  blade  and  = V/£R. 


Substituting  equation  B-2  into  B-l  and  noting  that  dr  = Rdx, 
the  expression  for  thrust  becomes 


- . b hRc  ( il)‘ 

T 3 


2 

J ( x ♦ ti  s i n - ; CL 


dx 


(B-3) 


Expressing  equation  B-3  in  terms  of  coefficient  of  thrust,  we 
have : 

A 


Cj- 


r-A  ( *RJ 


_ bcR 

7 ' "7T 


/ 


(x  + i*  sin  • ) dx 
x 


(B-4 ) 


From  equation  B-4  we  are  ready  to  calculate  the  thrust 
coefficient  as  follows: 


2bcR  .1 


= -r-  = ~rar/  (x  + li  sin  *)2cl  dx 

—7 — X 


■/ 


(x  ♦ - sin  . ) CL  dx 
x 


(B-5 ) 


Then  to  calculate  t„  for  each  u we  would  use  the  maximum 

cmax 

lift  coefficient  at  each  blade  station.  This  can  be  expressed 
as  follows:  ^ 


max 


/ 


(x  + p sin  •)  C,  dx 
tnax 

x 


(B-6 ) 
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From  Figure  21  the  retreating  side  of  the  blade  is  partially  in 
the  reversed  flow  region.  The  locus  of  this  region  is  defined 
by 


x = sin.  ( 180°  S * £ 360°) 

Then  equation  B-6  would  be  integrated  from  x = --  sin  . 
to  x = 1 . 

The  effect  of  pitch  rate  on  maximum  rotor  thrust  is  determined 
as  follows: 


M = 21^‘q  sin  > 


or  nondimensiona 1 ly  as 


M 


M _ 2q  sin 


< Xb 

The  total  rotor  moment  of  the  advancing  blade  if  expressed  as 
follows : 


(B-8) 


OTAL 


(x  + ,1  sin  •)  C 


. x dxl 
max  / 


+ - ( B-9 ) 


C uQ  R • 

where  V = — s , I,  is  the  blade  inertia  and  a is  the  lift 

d D 


curve  slope. 
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SIMPLIFIED  MAXIMUM  THRUST  PROGRAM 

/•  TC-I4U  UY  NUMERICAL  INTEGRATION  ANO  TABU  LOOK-UP  FOR  Cl-f(MACH)  • /; 

.PICUM  . iMXl  « LUi  « HfflOli  i 

aH2),a1U)a»; 

aMl),afitl)Mj 

«lp»|a.m>W92Sl,  

GET  LIST(OPT); 

IF  OPT-U  THEN  00  TO  »2; 

. GET  ILSKcinL 

GET  LI  ST(nunaa#|>r#R,om«f«r,C#  rho,v»,«,  Ihlarla); 

GET  LISTUK); 

PUT  LIST(");  , .... 

IF  KK*1  THEN  PUT  LISTE'FULL  UPLOAD  IN  REVERSE  FLOW  REGION' ) j 
IF  KK*-1  THEN  PUT  LI  T(  'FULL  DOWNLOAD  IN  REVERSE  FLOW  REGION' >j 

. If  AfU  THEM  PUT  Lis,  f MO  LOAU  Ik  mm  PLOW  AISIOM1)) 

PUT  L I ST( ' * ) ; 

PUT  LIST!  1 THE  FOLLOW  I NO  IS  A LIST  Of  THt  CL  DATA  TASLE  USED  FOK  COMPUTATION')) 

PUT  LIST!  'MACH  NO*  . OilAJL'lj 

OU  t«il  TU  2 II; 

PUT  IIIAnrU«.US,cln(Z*l))(  I HI  ■ 

-LUAGXi  . - 

END  Hi 

PUT  l ISTl  * * J;  

PUT  IIST("); 

I NU  • . ; 


a I (1),aM2),a1<  3), ahlN ) amtpj 
l)U  13-1  TU  It; 

I m-U;  

IF  MlaUAI  3»2|  1 3*N  THEN  GO  TO  ilk; 
a I I ■*  I M3); 

ihUiAliUi); 

ahlt-al  l ) * 

UU  U-l  TO  11; 

a««II*Ji*(U>11;  

m1»1; 

IF  I 3 >2  THEN  q 1 • a ; 

IF  13-2 lUik  THEM  aliflUU) 

nacl'»*iD»(  '|2  )*MOK; 

CALL  TAULElmadvcLclaJi 

I nt*lnt*cl •d*«ql*q2*q2; 

END  *6; 

(41  i ) ■ I nC; 

END  »!>; 

I ntb  *0; 

1 n (b*U;  — 

l3ba|(3)al(b)»K; 

DO  lb-1  TO  11; 

*■•!*( 1G-1 ); 
ql  up; 
i iacii*M()U*ql; 

CALL  TABLE 

q2  l*d  •>ll*ql; 
lntS"lnt!»+«|2; 

— 4M»alMb»«2a<*»  — 

IF  I ntb>  I 3 4 THEN  GO  TO  »10; 

END  &9; 

lnt$a|ntb-<|2»  

Intb»lntb-q2#*; 
q3»l 3U-lnt6; 


xolri-xU; 

A-*O*«|3/(q2*A-l0); 

IF  «-a0*U  THEN  GO  TO  

t|l»«*n*p; 

t-tachaMOK*ql ; 

CAU  TABLE(t»aaNy#l,<(ai>» 

i|2a(ft-xO)ae1«Ml«<tl) 
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i. 

i. 


L 


INPUT  COEFFICIENT  PROGRAM 


/•  INPUT  COEFFICIENT  DETERMINATION  FOR  AFSA  CONTRACT  •/; 

DECLARE  hp'.m  , hpl(7)  , hptl(7)  , hpll<7)  , vc(7)  , hpcl(7)  # hpcl<7>  j 

GET  LIST(n); 

GET  LI$T(hpl,hp2, vc)} 

GIT  

GET  L I ST( h,o*t ) ; 

GET  LIST(twl,s*2)j 



t«o«t*273. If; 
th*t* "t/211 .11; 



tl<p*(ltt4/thPtaj 
rho«s I £p*« 002)71; 

iWtWNWi 

•"6.2*; 

•r*p*pl • p**2 ; 

§•  t *<*•«/(  pi  •#  >t — 

nondf  •rho*prtp«wr**2j 
rwndp«rx>ndf  •wr/550; 

■ 

ctl«|wl/nondf j 
ct2»(jw2/nondf ; 

/•  CALCULATION  Of- INOWOtO  PR*|A  fl 

cti»*ctl; 

•ll  00  l«l  TO  m; 

- I 

mu "v/wr; 

CALL  vlnd; 



IF  mu). 14  THEN  ki»lj 
hplll  l)*khctp»nondf*»l/550j 

■ gNO  «1; 

ctp»ct2; 

•2 1 DO  I "l  TO  mj 



mu«v/wr; 

CALL  vlnd; 

fHWl 

IF  mu > . 1 4 THEN  kl»l* 

hpl  2(  I )«kl*cti»*r>ondf  "v 1/550  j 

€N0  el ; 

•)t  DO  I«1  TO  m; 

hpcl( I )"hpl(  I )-hpl 1 ( I ) j 

hpcM+HNpH+»  hplHI»; 

END  «); 

/•  CALCULATION  UF  DELTA  2 •/; 

t-lj  

dhp"hpc2 ( I ) -hpc 1 ( I ) ; 
v«vc( I )"1.6*7*j 

■ uiu"vfwrj  1 " 

ixu2amu*tHj; 

kl"«  lg"nondi>"(  l*4.6"mo2 )•(  ctc»ctl )"«2/l; 

k2»*lu*nondp*<  ; 

d«  1 2 ><tli|)/  (k2*kl ) ; 

/•  CALCULATION  UF  HP  STALL  •/; 

Twrr 

v"vc(r»)*1.6*7*; 

wu*v/wr; 

mu2*nu»iHi;  * 

hpd»d«l2*»lc*f>ondp»(X*4.G*mu2)«(etl*ctc)**2/l; 

hpcl(m)«hpcl<m)-hpd; 

hpd-del  2MlB*nondp*Cl*4.6*mu2)*Cct2*cte)««2/0j 

hpc2  (i  0 »hpc2  Cm)  -Hpd ; 

tci-2*cti/»ic; 

tc2»2*ct2/*Tfi;  

leer  t ■.  l*.2/s<irt(  1*5  Q*mu2); 

IF  tcKtccrt  THEN  GO  TO  «k; 

dqttr."tcl,tecrt ; ' 

lip«"200U0*d«ltc**l.S; 

hpcl(m)*lipcl(n)-hpi; 

oo  to  *5; 

•4:  hpt"U; 

•St  IF  ec2<tccrt  THEN  GO  TO  ••; 

0iire»cc2-tcertj  ■ - 1 1 ' 

hp»"2UUUU"dtl  tc**l.S; 
hpc2(n)*hpc2(M)-l»p»; 

GO  TO  *7  ; 

•it  hp*»U; 

/*  CALCULATION  UF  mcro  •/; 

rn np«nnc2(nj*n]icrtnrT 

IF  hp>U  THEN  fit)  TO  til) 
mcro*,/  ■ 
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00  TO  *2A; 

•2)1  dm*. 7S«(tc2-tcl)j 

di>«hn*  100000/ nondpj 


k*0; 

mcroa.7N;  

§Ti  ” mcromncro-,  005J 
•St  d«lni»Untu)«wf/¥»-mcro*.  7S*tcl) 
d«lm2a(l*mu)»wr/v»incro*«  7S«tc2; 


dcp2ad«1n2**3*l  . UQI)-d«1m2»(  . Q22-.ll*d«1m2)} j 
dcpladcpl*10UQQQ; 

dep2adCp2»iTOO0Oj  

dcpdn>(dcp2-dcpl )/  (d*lm2-dtlml ) ; 

IF  kal  THEN  00  TO  •*; 


if  nnawMcpawro  twwr  oo  Tint) 

OO  TU  §2; 

• ll  dcpdnl adcpdm; 

uKromicroi  .005  J 

k-lj 

OO  TU  I); 

-in wyFgTTJenJn*qDaw?w.oo?rcuLpoiwfUuP>ii)iiA.iuj 

dcpl*dcpl*nondp/lUQUUU j 

hpel(n)  ■hpel<M)*dcpl) 

/•  CALCULATION  OF  d*10  *t\ 

• 2 A t v *vc ( 2 )• 1 • OH  71 ; 

MHi*w/«r  j 
»3-(X*nu)»wfAvi] 

IF  a J<»*cro  THEN  00  TO  «21j 
de ln»*J*ncro*, 75*(2*ct2/tlg); 
dcp-.l«li  !••)•(  ,0033-dtTn*C  .022-,  HMdtPT)  7) 
dhpadcp*nondp; 

hpcl(2)ahpcl(2)*dhp; 


“SIT  1 mu2*i'lu*Mu; 

•li  hpf  ■45*rho*v*v*v*fo/SS0j 

' hPc.liocl(2J"l»pf“NpdlJ  * 

d.lU-hpe/kl; 


/*  cmoprmm  ut  r »rr 

VVC<m)«l.M7l; 

imi'v/wr; 

mu2"inu*mu; 

hpd-d.l0Ml(*nondp«< 

hpfe-lipcllml-hpd; 

kl".5*rho*v«v«»/S50| 

f.hpfc/kl; 

_/•  CALCULATION  OF  TOTAL  POWIA  «/> 
k-0; 
et*ctl; 

00  I«1  TO  m; 


•111 

jUx_ 


• »cTTT>nn7Tj 

mu»v/wr; 

etp*«uH(«A*i.k/(nondf»nondf  )»el»et>» 
lindl»d<:IO»>l»«nondp«(l»>.t»wul)/l 


hpd2*d«l  2 • » I .•oondp»(  1*1. 
hp«»»k«v/5SQ; 

tccrt>.l*.2/*»|Tt(  J7 

IF  tcUrert  THEN  00  TO  till 
d«l te»te-tccrt 


mf 


hp».2uo<)u*d«l  t e • • 1 . 5 ; 

SO  TO  .12; 

_hp»*U; 


• 111 


*}.<l»Mu)*><r/y»; 

IF  x3>mcro  THEN  GO  TO  ilk; 


GO  TU  .13; 

.IN  I d.lm.«3-mcro».75*tc; 

dca.d.ln«»>«(.00»-d.lw«(.022-.ll«d«lw>>) 

•15 1 hpc.'lep*nondp; 

CALL  vlnd; 
kl-l-(nU-.lk).k3 


IF  i iu> t Ik  THEN  kl«I; 
hi>3»kl*ctp*fu>ndf  •vl/550j 

hot* hndl*hnd2»hP3*hak*hp»*hpcj 

IF  k-1  THEN  GO  TO  •16; 


•rror*hpl(  U-hpt; 
tsi Cl 


air&*irrar*lHg/hflllJL 

PUT  IMAGC(<Ml.vc(l)»h| 


,hpl(  t ),hpt,*  ror,p«rcH TWTT7 
trrof » 


PUT  IMAGE(fwlrvc( 
I ml | IMAGE; 

— » _Jt» • — JlflV 

GO  TO  «17; 
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<fl>« 


• rror «hp2( O-hPt  j 

a»fc-«ffflt«iia/hilll ji 
n>T  iMAflKtwJ, veil), hen 
IHAOE J 

ho»  h» 

(NU  ill 

IF  k*l  THIN  00  TO  *10/ 


dill*  


*k«.}»rho«**v*foj 
ml •*»•*/ (e»p*nondf)j 

xl'Ctp*  nondfV  (1*  fha«Af  u } j 

IF  THIN  00  TO  »l7| 

vl»1.01l*»«rtl«2)j 
IF  vbQ  THEN  00  TO  Mill 
*Io»mI/vj 

»l •Ml«Ml/(v*»*,IO»*»lH«*ln)**. Jj 

j*  ,H' thin  no  TO  Mill 

«ount*eount«li 

IF  counOkO  Then  00  TO  llOi 

mIm.-«mmI»-1mmIiii 

put^iist! 'calculation  non  convinoent  induced  velocitt  to  o')i 

nItOi  

00  TO  Mill 
vl««l*lll 

II  ..1.1ATA  THAN  ..«■ 

(NO  Vlnd) 
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